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Abstract. We develop new solvability methods for divergence form second order, 
real and complex, elliptic systems above Lipschitz graphs, with L 2 boundary data. 
The coefficients A may depend on all variables, but are assumed to be close to 
coefficients A$ that arc independent of the coordinate transversal to the boundary, 
in the Carleson sense \\A — Aq\\c defined by Dahlberg. We obtain a number of a 
priori estimates and boundary behaviour results under fmiteness of \\A — Aq\\c- 
Our methods yield full characterization of weak solutions, whose gradients have L2 
estimates of a non-tangential maximal function or of the square function, via an 
integral representation acting on the conormal gradient, with a singular operator- 
valued kernel. Also, the non-tangential maximal function of a weak solution is 
controlled in L2 by the square function of its gradient. This estimate is new for 
systems in such generality, and even for real non-symmetric equations in dimension 
3 or higher. The existence of a proof a priori to well-poscdness, is also a new fact. 

As corollaries, we obtain well-poscdness of the Dirichlct, Neumann and Dirichlct 
regularity problems under smallncss of \\A — Aq\\c and wcll-posedness for Aq, 
improving earlier results for real symmetric equations. Our methods build on an 
algebraic reduction to a first order system first made for coefficients A by the two 
authors and A. Mcintosh in order to use functional calculus related to the Kato 
conjecture solution, and the main analytic tool for coefficients A is an operational 
calculus to prove weighted maximal regularity estimates. 



MSC classes: 35J55, 35J25, 42B25, 47N20 

Keywords: elliptic systems, maximal regularity, Dirichlet and Neumann problems, 
square function, non-tangential maximal function, Carleson measure, functional and 
operational calculus 

1. Introduction 

In this article, we present and develop new representations and new solvability 
methods for boundary value problems (BVPs) for divergence form second order, 
real and complex, elliptic systems. We look here at BVPs in domains Lipschitz 
diffeomorphic to the upper half space R+ +n := {(t, x) G R x R n ; t > 0}, n > 1. 
The same problems on bounded domains Lipschitz diffeomorphic to the unit ball, 
contain noticeable differences, such as use of Fredholm theory, which we address in a 
forthcoming paper [6]. Here, we focus on the fundamental scale-invariant estimates. 

The system of equations is 



n rn 

(1) Lu°(t, x) = Y,Y, d * { A Zf& *) V(*> *)) = 0, 

i,j=0 13=1 

in R+ +n , where do = and di — ^7, 1 < i < n. We assume 

(2) a = (Aif(t, ./•));:/ e ^(R 1+n ; c(c^n), 



a = 1, . . . , m 
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and that A is strictly accretive on TL, meaning that there exists k > such that 

n m „ n m „ 

( 3 ) EE MAtf(t,x)f^x)Jt(x))dx>Kj2J2 \f?(x)\ 2 dx, 

for all f en and a.e. t > 0. The definition of ft, a subspace of L 2 (R n ; C^ 1+ ^ m ), 
will be given in Section |2j 

We seek to prove well-posedness for ([TJ, i.e. unique solvability in appropriate 
spaces given Dirichlet data u\ t= o, Neumann data d UA u\ t= o or Dirichlet regularity 
data V x u\ t= o, assumed to satisfy an L 2 condition. Note that the continuity estimate 
required for well-posedness in the sense of Hadamard is not included in our notion of 
well-posedness, but will be shown to hold. For the Neumann and Dirichlet regularity 
problems, we will work in the class of weak solutions whose gradient Vt, x u has L 2 
modified non-tangential maximal function N*(Vt, x u ) m L 2 . (See Definition 12. 11 ) 
Under our assumptions, we shall describe the limiting behaviour of Vt, x u at t = 
and oo and obtain a perturbation result for well-posedness in this class. For the 
Dirichlet problem, it is more natural given our method to work in the class of 
weak solutions with square function estimate JJpi+n | Vt, x u\ 2 tdtdx < oo. Under our 

assumptions, we shall describe the limiting behaviour of u at t = and oo and prove 
a rigidity theorem that shows new a priori non-tangential maximal estimates and 
L 2 estimates, and obtain a perturbation result for well-posedness. 

Let us begin by pointing out that the coefficients depend on t, which makes 
these problems not always solvable in such generality. In Caffarelli, Fabes and 
Kenig [T5], the necessity of a square Dini condition is pointed out. There has been 
a wealth of results for real symmetric equations (i.e. m — 1 and A+j = Aji G R, 
% = L 2 (R n ; C 1+n )). In Fabes, Jerison and Kenig [26], the L 2 Dirichlet problem is 
solved under the square Dini condition and continuity. Dahlberg removed continuity 
and proved in [19] that if the discrepancy A\ — A 2 of two matrices A±, A 2 satisfies a 
small Carleson condition, i. e. if \\A\ —A 2 \\c from Definition 12. II is small enough, then 
/^-solvability of the Dirichlet problem with coefficients A-y implies L p2 -solvability 
of the Dirichlet problem with coefficients A 2 with p 2 — p\. The smallness condition 
was removed in Fefferman, Kenig and Pipher [28], but then the value of p 2 becomes 
unspecified. R. Fefferman obtained in [27] the same conclusions as Dahlberg with 
p 2 — pi, under large perturbation conditions of different nature. See also Lim pET] . 
Kenig and Pipher [33] proved that the L p -Neumann and regularity problems are 
uniquely solvable if the discrepancy A(t, x) — A(0, x) satisfies Dahlberg's small Car- 
leson condition, depending on p G (1,2 + e). Moreover, in [34] they proved small 
perturbation results for the Neumann and regularity problems analogous the result 
[T9] for the Dirichlet problem, as well as large perturbation results for the regularity 
problem analogous to [2H] for the Dirichlet problem. 

Some related results of Kenig and Pipher [33] (going back to questions of Dahlberg 
|19j). Dindos, Petermichl and Pipher [21] and Dindos and Rule [2Sj are obtained 
under smallness of a Carleson condition on tVt, x A(t, x). See also Rios' work [4*2] . 
Such an hypothesis does not compare to the one on A(t,x) — A(0, x). 

We note that these results are obtained for L p data, for appropriate p's, including 
p = 2. This is using all the available technology for real scalar equations, starting 
from the maximum principle, hence L-harmonic measure, and Green's functions. 
Moreover, as far as solvability is concerned, the main thrust of these works is to get 
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p = 2 with non-tangential maximal estimates, using for this the classical variational 
solutions, or those obtained via the maximum principle. 

Of course, t-dependent coefficients incorporate the t-independent ones. We refer to 
the book by Kenig [30] and references therein, and to Alfonseca, Auscher, Axelsson, 
Hofmann and Kim [I] for more recent results on perturbation of real symmetric 
(or complex and constant) equations. See also below. 

We mention a series of works for two dimensional equations on the upper half- 
plane with t-independent coefficients. Auscher and Tchamitchian [11] study complex 
coefficients equations with diagonal A (which we call here block form) and describe 
Dirichlet, regularity and Neumann problems for L p data for p > 1 and even for 
data of Hardy type for p < 1. This is a precursor of the work for systems here, 
as it built upon new proofs relying on Calderon-Zygmund operators (which are 
no longer available here) of the one dimensional Kato conjecture proved earlier by 
Coifman, Mcintosh and Meyer [16] and its extension by Kenig and Meyer [32]. For 
real equations of non block forms, Kenig, Koch, Pipher and Toro [31] show that the 
Dirichlet problem is well-posed for large enough p (and obtain counterexamples for 
any given specific p) by showing that L-harmonic measure is absolutely continuous 
with respect to Lebesgue measure on the boundarj^. Kenig and Rule [SB] then 
obtained well-posedness for the Neumann and regularity problems with p — 1 > 
small enough (and obtained counterexamples for any specific p > 1). The recent 
thesis of Barton [H] deals with complex, t-independent perturbations of the 
situation in [31] [36] . and she obtains well-posedness of the Neumann and regularity 
problems in L p for p — 1 > small and even at p = 1 with data in the classical 
Hardy space. 

As the reader has observed, we consider complex systems and we wish to obtain 
L 2 solvability under conditions as general as possible (we mention that L p solvability 
with our methods when p ^ 2 is under study at this time). For this, we need other 
tools than those mentioned above. In fact, the tools we develop and that we describe 
next would not have been conceivable prior to the full solution in all dimensions of 
the Kato problem and its extensions. In Auscher, Axelsson and Mcintosh [9] , a new 
method was presented for solving BVPs with t-independent coefficients, following an 
earlier setup designed in Auscher, Axelsson and Hofmann [7]. The main discovery 
in [9] is that the equation ([Q) becomes particularly simple when solving for the 
conormal gradient defined by 



(4) / = V A u 



d VA u 



where d VA u denotes the (inward for convenience) conormal derivative (see Section HI), 
instead of the potential u itself. It is a set of generalized Cauchy-Riemann equations 
expressed as an autonomous first order system 

(5) d t f + DBf = 0, 

where D is a self-adjoint (but not positive) first order differential operator with 
constant coefficients that is elliptic in some sense and B is multiplication with a 
bounded matrix B(x), which is strictly accretive on the space % in (jH]) and related 
to A(x) = A(t,x),t > 0, by an explicit algebraic formula. The operator DB is a 



1 A recent preprint by Dindos, Kenig and Pipher [23] posted during the revision of this article 
shows this is related to well-posedness with BMO data. 
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bisectorial operator and can be shown to have an L 2 -bounded holomorphic func- 
tional calculus for any (t-independent) matrix A satisfying (j2J) and (jBJ. This fact 
was proved earlier by Axelsson, Keith and Mcintosh [T3"l Theorem 3.1] elaborating 
on the technology for the solution of the Kato problem by Auscher, Hofmann, Lacey, 
Mcintosh and Tchamitchian [10]; a more direct proof is proposed in Auscher, Axels- 
son and Mcintosh [8]. As explained there, the main difficulty is the non-injectivity of 
D. The upshot is the possibility of solving (JSJ) by a semi-group formula / = e~*' DS '/o 
with /o in a suitable trace space, and such / has non-tangential and square function 
estimates. The BVP can then be solved in an appropriate class if and only if the 
map from the trace functions to boundary data is invertible. This is the scheme 
for the Neumann and regularity problems, for which the boundary data is simply 
the normal or tangential part of VaU- For the Dirichlet problem, it turns out that 
a "dual" scheme involving the operator BD can be used similarly. The one-to-one 
correspondence between trace functions fo and boundary data may fail, see Axels- 
son |12| . and it is here that restrictions on A appear. It is known to hold if A is 
(complex) self-adjoint or block form (i.e. no cross derivatives d^Adi or diAd$, i>l, 
in (CQ)), or constant. Another consequence of this method, and this is why consid- 
ering complex coefficients is useful, is that the set of t-independent A's for which 
solvability holds is open in L^. See [9]. 

Our work for t-dependent coefficients takes the algebraic reduction to (GO) as a 
starting point, the conormal gradient becoming the central object. We shall state 
the main results in Section [2] and explain the strategy in Section [3J It involves in 
particular study of a highly singular integral operator Sa, with an operator- valued 
kernel. On a technical level, proper definition and handling of this operator is most 
efficiently done using operational calculus rather than the usual maximal regularity 
treatment originally due to de Simon [22] (see Kunstmann and Weis [38l Chapter 
1] for an overview) and this avoids having to assume qualitatively that A is smooth 
in the calculations. We use the terminology operational calculus, following the 
thesis [1] of Albrecht, for the extension of functional calculus when not only scalar 
holomorphic functions are applied to the underlying operator (in our case DB with 
Bq(x) = B(0,x)), but more general operator-valued holomorphic functions. The 
Hilbert space theory we use here to prove boundedness on appropriate functional 
spaces in Section [7J surveyed in Section 16. 1\ is a special case of the general theory 
developed in Albrecht, Franks and Mcintosh [Bl Section 4], Lancien, Lancien and 
LeMerdy [39] , and Lancien and LeMerdy [ID] . For further details and references, we 
refer to Kunstmann and Weis [3E1 Chapter 12]. 

The Carleson control on the discrepancy A(t,x) — A(0,x) from [TBI EH [331 [M] 
appears in a very natural way in the estimates of Sa, and well-posedness of the three 
BVPs with coefficients A(t,x) will follow under smallness of this control and well- 
posedness for coefficients A(0,x). We mention that the Dirichlet problem could be 
obtained by an abstract duality procedure from a regularity problem, in agreement 
with the results of [331 [31] f° r rea l symmetric equations. See also Kilty and Shen [37], 
and Shen [IB]. We will formalize this abstract procedure in our subsequent work [6]. 
We remark however that although the hypotheses are the same for each BVP, the 
perturbation results can be proved independently of one another. For example, one 
does not need knowledge on well-posedness of regularity for A or of Dirichlet for A* 
to obtain well-posedness of Neumann for A, in contrast with the results in [34"] . 
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We do not know how to prove well-posedness under the finiteness of \\A(t, x) — 
A(0, x) ||c only. (In light of [2E1 EH], this would first require to extend our methods to 
solvability for L p data.) However, thanks to our representations, we do obtain under 
this hypothesis a number of a priori estimates and boundary behaviour on solutions 
to the equation (CQ) without any assumption on well-posedness. For example, and we 
concentrate on this to finish this introduction, we show that if \\A(t, x) — A(0, x)\\c < 
oo, all weak solutions to (CQ) with coefficients A such that // R i+« | Vt, x u\ 2 tdtdx < 

oo satisfy, up to a constant, sup-L2 estimates as well as non-tangential maximal 
functions estimates. More precisely, we obtain 

max(sup ||u t || 2 , IW^Ih) < || Vt,xU\\L 2 (tdt;L 2 ) 
t>0 

provided the right hand side is finite and the solution u vanishes at infinity in some 
sense (see Section [5] for precise meaning). Note in particular that this applies when 
A is t- independent and in that case, this is implicit from [HI Corollary 4.2] when re- 
stricted to the class of solutions considered there. Domination of the non-tangential 
maximal function ||iV*(w)||2 by the square function \\V t,xu\\L 2 {tdt-,L 2 ) ~ II^MIh) 
S(u)(x) = (J| y _ x | <f |V 't t yu\ 2 dtdy /i™ -1 ) 1 / 2 , is reminiscent of the result of Dahlberg, 
Jerison and Kenig [20], and also of Dahlberg, Kenig, Pipher and Verchota [21]. But 
there is a difference. In [20] comparability of N*(u) and S(u) is obtained for solu- 
tions of the equation ([T]) under ([2]) and ©, A real and m — 1, in all L q (R n ; dfi) 
spaces, < q < oo, with \i a doubling A ra weight with respect to L-harmonic 
measure. If the Dirichlet problem in the class ||A^,(u)|| p < oo is proved to be solv- 
able for one 1 < p < oo, then Lebesgue measure is A^ of L-harmonic measure, 
hence ||iV*(ii)|| g ~ || >!?(«) ||,j. This fact follows in particular from combining |29j and 
[2Sj under \\A(t, x) — A(0,x)||c; < oo and A real symmetric. In [21], comparability 
||iV*(ti)|| g ~ < q < oo, is obtained for real symmetric constant elliptic (in 

the sense of Legendre-Hadamard) second order systems (and also higher order but 
the formulation becomes different) on bounded Lipschitz domains owing to the fact 
that L 2 solvability of the Dirichlet problem was known (see the introduction of [2~T]). 
This comparability also follows for real non-symmetric scalar equations in two di- 
mensions combining the results of [HI] and again [20]. Here, although we obtain only 
one part of the comparison, it is essential to note that this is an a priori estimate 
valid independently of well-posedness. The existence of an a priori proof is new even 
for real symmetric scalar equations under the Carleson control (for example for all 
the t-independent ones), and is permitted by the solution of the Kato square root 
problem and its extensions. 
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2. Statement of results 



In this section we state our results concerning solvability of boundary value prob- 
lems on the half space R 1+n , and show how they extend to domains Q which are 
Lipschitz diffeomorphic to R + +n . 

Let us first fix notation. We write {eo, e±, . . . , e n } for the standard basis for R 1+n 
with e "upward" pointing into R+ +n , and write t = x for the vertical coordi- 
nate. For the vertical derivative, we write do = d t . For an m-tuple of vectors 
v = ( t, j Q )o<i<^n ?1 ' we wr he v x and for the normal and tangential parts of v, 
i.e. (fj_)o = anc ^ i v ±)? = when 1 < % < n, whereas (f||)f = when 1 < i < n 
and (f||)o = 0. We write ft(x) := f(t,x) for functions in R^_ +n . As compared to 0, 
we here use subscript to denote restriction to the boundary R n at t = 0, rather 
than the normal component of /. We also prefer to use small letters f,g,... to 
denote functions in R+ +n , since this is where we work most of the time, not on the 
boundary as in [9]. 

For tuples of functions and vector fields, gradient and divergence act as (Vt. x u)" = 
diU a and (div t ,xf) a = Y17=o^ifi '■> with corresponding tangential versions V x u = 
(Vt i3; M)|| and (div x f) a = Y^H=i®if?- With curl^/ = we understand that djf? = 
dif®, for all i, j — 0, . . . ,n. Similarly, write cuY\ x f n = if djf° = dif° ! , for all i, 
j = 1, . . . ,n. 

We assume that A is strictly accretive on H, i.e. 



holds for some constant k > 0, uniformly for t > and all / belonging to the closed 
subspace 



For systems, ([3]) is stronger than a strict Garding inequality on R 1 ^™ (i. e. integration 
would be on R+ +Tl and / such that curl^/ = 0); still ([3]) is natural given the type 
of perturbation we consider here. 

The boundary value problems we consider are to find u G P'(R^ +n ; C m ) solving 
the divergence form second order elliptic system 



in distributional sense, with appropriate interior estimates of V t>x u and satisfying 
one of the following three natural boundary conditions. 

• The Dirichlet condition u = p on R n , given p> e L2(R n ; C m ). 

• The Dirichlet regularity condition V x u = <p on R n , given p e L2(R n ; C nm ) 
satisfying curl^ = 0. 

• The Neumann condition (eo, AVt, x u) — tp on R n , given p e L 2 (R n ; C m ). 

Note that since we shall impose distributional Vt^w G L l £ c , u can be identified with 
a function in W 2 1,loc , i.e. with a weak solution. 





<Hv ttX AV t ,x u — in R 
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Definition 2.1. The modified non-tangential maximal function of a function / in 
R^_ +n is 

N.(f)(x) := sn P t-^ 2 \\f\\ L2mttX)h x G R n , 
t>o 

where W(t, x) := (cq H, cot) x B(x; cit), for some fixed constants cq > 1, ci > 0. The 

1/2 



modified Carleson norm of a function g in R+ +n is 



Iff , , 2 dtcfe 



MIc := | SU P7777 / / SUp |0, 

3 M JJ(0,l(Q))xQ W(t,x) 1 

where the supremum is taken over all cubes Q in R n , with l(Q) denoting their side 
lengths. 

Note that different choices for cq, c\ will give different, but equivalent norms 
ll-^*(/)l|2) as weu as equivalent norms \\g\\c- Furthermore, this maximal func- 
tion is really non-tangential since iV*(/) and the closely related maximal function 
sup|j,_ x i <t t~^ 1+n ^ 2 \\f\\L 2 (w(t,y)) have equivalent L 2 norms. The latter was introduced 
in [33] • The modified Carleson norm originates from Dahlberg [TP] . 

We will use the modified Carleson norm to measure the size of perturbations 
of t- independent coefficients Aq. (In fact we shall use a possibly weaker quantity 
|| • ||*; see Section 0) Intuitively, \\A — Aq\\c < oo means that in a certain sense 
A(t, x) = Aq(x) at t — 0, but also that A(t, x) is close to Aq(x) at all scales since we 
are dealing with a scale invariant norm. Also, given A, such a Aq is unique and has 
controlled bounds thanks to the following lemma. We remark here that the modified 
Carleson norm there can be weakened to the usual one. See Section for proofs. 

Lemma 2.2. Fix A : R^_ +n ->■ C(C^ +n ^ m ) with < oo and strictly accretive on 

%, with constant of accretivity k > 0. Assume that Aq are t-independent measurable 
coefficients such that \\A — Aq\\c < oo. Then Aq is uniquely determined by A, i.e. 
if A' are t-independent coefficients such that \\A — A' Q \\c < oo, then A' Q = A almost 
everywhere. Furthermore A is bounded and strictly accretive on H, with 

« < ^0 < Po||oo < H^IU, 

where k,q denotes the constant of accretivity of Aq. 

For the Neumann and Dirichlet regularity problems, our result is the following. 

Theorem 2.3. Consider coefficients A G L 00 (R^ +n ; £(C^ 1+ri ^ m )) which are strictly 
accretive on TL. 

(i) A priori estimates: Assume that u G VK 2 1 ' Ioc (R^ +n , C m ) has gradient Vi jX w 
with estimate || ( V t || 2 < oo, and where u satisfies (QP in R^ +n distri- 
butional sense. If there exists t-independent measurable coefficients A such 
that \\A — Aq\\c < oo, then V^w has limits 

limir 1 \ \\V s>x u s - g \\lds = = lim t~ l / || V sx u s \\lds, 

for some function go G /^(R"; C ( - 1+n ^ m ), with estimate 1 1 t?o 1 1 2 ^ || iV* ( V t ;r w) || 2 . 

(ii) Well-posedness: By the Neumann problem with coefficients A (or Aq) being 
well-posed, we mean that given tp G L2(R n ; C m ), there is a function u G 

unique modulo constants, solving (QJ), with coefficients A 
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(or A ), and having estimates as in (i) and trace go = lim^o Vt )X iz such that 
(A g ) ± = (p. 

The following perturbation result holds. If the Neumann problem for A 
is well-posed, then there exists e > such that if \\A — Aq\\c < e, then the 
Neumann problem is well-posed for A. 

The corresponding result holds when the Neumann problem is replaced by 
the regularity problem and the boundary condition (Aogo)± = <p is replaced by 
(<?o)n = ¥ £ L 2 (R n ; C nm ), where <p satisfies curL^ = 0. Moreover, for both 
BVPs the solutions u have estimates 



(iii) Further regularity: Assume that A is as in (i), with \\A — A \\ c sufficiently 
small and consider solutions u as in (i). 

If A satisfy the t-regularity condition \\tdtA\\c < oo, then 



t H> Vt, x Ut £ L 2 is continuous and lim t _> II V ttX u t -g \\ 2 = = limbec || V t , x u t \\ 2 . 
The converse estimate \\N^(Vt, x u)\\l < / °° ||o^Vt, x w|||tdt holds provided \\td t A\\c 
is sufficiently small. 

If max(\\tdiA\\c, \\td t A\\c) < oo holds for some i — 1, . . . , n, then 



The estimate \\N„{y t , x u)\\\ < Yli=i Jo°° ll^i^t^^lli^ holds provided \\£Vt, x A\\ c 
is sufficiently small. 

Under the hypothesis \\A — Aq\\c < oo, the implicit constant in (i) depends on n, 
m, 1 1 ^4 1 1 oo , k. In (ii) and (iii), under the smallness hypothesis on ||v4 — AqIIcs which 
depends on n, m, WAW^, n (and in (ii) also on the "well-posedness constant" for A ), 
the implicit constants depend on n, m, ||v4||oo ; K (and in (iii) also on the regularity 
assumptions on A). However, some inequalities are true without smallness and/or 
well-posedness as the reader can check on reading details in Section [HJ Note that in 
(ii), the uniqueness holds in the class defined by ||A^(Vt )a; tt)||2 < oo. 

For the Dirichlet problem, our main result is the following, including a rigidity 
result of solutions with square function estimates. 

Theorem 2.4. Consider coefficients A £ L 00 (R l ^' n ; £(C^ 1+n - )m )) which are strictly 
accretive on H. 

(i) A priori representation and estimates: Assume that u £ W 2 ' oc (R+ +n , C m ) 
has gradient bounds J °° || Vt, x u t \\ 2 tdt < oo and satisfies (T7]j in R^ +n distri- 
butional sense. If there exists t-independent measurable coefficients A such 
that \\A — Aq\\c < oo, then u = u + c almost everywhere, for a unique 
u £ C(R + ;L 2 (R n ; C m )) and constant c £ C m . Identifying the functions u 
and u + c, it has L 2 limits 



N*(Vt,xu)h ~ 1Mb « 1Mb- 




d t V t , x u\\ltdt < sup || V t)a! ut Hi « ll-^*(Vt )a! u)||l, 




lim \\ut — Mo — c| 1 2 = = lim \\u t — c\\ 2 , 
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for some u £ L 2 (R n ; C m ), and we have estimates 

POO 

max( 1 1 N*(u) \ \ 2 , sup 1 1 ut || 2) < / || V t , x u\\ltdt. 
t>o Jo 

(ii) Well-posedness: By the Dirichlet problem with coefficients A (or Aq) being 
well-posed, we mean that given ip £ L 2 (R n ; C m ), there is a unique u £ 
W 2 oc (R+ +n , C m ) having gradient bounds and solving ([!]), with coefficients A 
(or A ), as in (i), and trace uo = <p. 

The following perturbation result holds. If the Dirichlet problem for Aq 
is well-posed, then there exists e > such that if \\A — Aq\\c < e, then 
the Dirichlet problem is well-posed for A. Moreover, these solutions u have 
estimates 

POO 

||iV*(u)||i » sup \\u t \\l « / ||V M it||2^ « \\(p\\l. 

t>0 Jo 

Say that a function w £ L 2 oc (R + , L 2 ) vanishes at 00 in L 2 sense if linu._5.00 1 1 tt^t 1 1 2 = 

0. For the solutions u as in (i), we see that the following three statements are 
equivalent: u vanishes at 00 in L 2 sense, u £ C(R + , L 2 (R n ; C m )), u £ L 2 (R n ; C m ). 
So in (ii), the boundary condition uo £ L 2 (R n ; C m ) forces c = and u = u. 

Under the hypothesis ||A — Ao\\c < 00, the implicit constants in (i) depend on n, 
m, H^Hoo, k. In (ii), under the smallness hypothesis on \\A — A \\c, which depends 
on n, m, \\A\loo, k and the "well-posedness constant" for A , the implicit constants 
depend on n, m, \\A\loa, k. However, some inequalities are true without smallness 
and/or well-posedness as the reader can check on reading details in Section |HJ Note 
that in (ii), uniqueness holds in the class defined by / °° ||Vt )X u|||^ < 00. 

As mentioned briefly in the introduction, the hypothesis on well-posedness with t- 
independent coefficients A is satisfied, for all three BVPs, for Hermitean coefficients, 

1. e. A (x)* = Aq(x), for block form coefficients, i.e. (Ao)j_\\ = = (Ao)\ ]± , and for 
constant coefficients, i.e. Aq(x) = Aq, as well as for sufficiently small t-independent 
L 00 (R n ;£(C (1+ri)m )) perturbations thereof. This was proved in [HI Theorem 2.2]. 
That the notions of well-posedness of these BVPs used in P coincide with the ones 
here, for ^independent coefficients, follows from Corollaries 18.41 and 19.41 

Note that we do not assume pointwise bounds on the solutions, hence we use iV* 
instead of the usual non-tangential maximal function. 

When m — 1 and A, A are real symmetric (and R^ +Tl replaced by the unit ball), 
Theorem I2.3( ii) is in [33], and Theorem 12.4( h) is in [TS] (and [20] for the square 
function estimate). The rest of Theorems 12.31 and 12.41 is mostly new. In Section [31 a 
more detailed road map to the proofs is given. 



Proof of Theorems \2.3\ and\2.4\ For the Neumann and regularity problems in R 



1+?) 



part (i) follows from Theorems 18.21 and 18.31 part (ii) follows from Corollary 18.61 and 
part (iii) is proved in Theorem 1 1U. 41 

For the Dirichlet problem in R+ +n , part (i) follows from Theorems 19.21 and 19.31 
and part (ii) follows from Corollary 19.51 except for the estimate of the non-tangential 
maximal function, which is proved in Theorem 110.11 □ 

We end this section with a remark on the Lipschitz invariance of the above results. 
Let f2 C R 1+n be a domain which is Lipschitz diffeomorphic to Ri +n , and let p : 
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R^™ — > f2 be the Lipschitz diffeomorphism. Denote the boundary by X := d£l and 
the restricted boundary Lipschitz diffeomorphism by po : R n S. 

Given a function u : f2 — >■ C m , we pull it back to m := « o p : R^ +n — >■ C m . 
By the chain rule, we have V t , x u = p*{St,xU), where the pullback of an m-tuple 
of vector fields /, is defined as p*(f)(x) a := p t (x)f a (p(x)), with p l denoting the 

transpose of Jacobian matrix p. If u satisfies divt^AVt^u = in Q, with coefficients 
A G £(C( 1+n > m )), then u will satisfy div^AV ' t>!B u = in R^ +n , where A G 

L 0O (Ri +n ;£(C< 1+n > m )) is defined as 

(7) A(x) := |J(p)(x)|(p(x))" 1 i(p(x))(p*(x))- 1 , x G R^. 

Here J(p) is the Jacobian determinant of p. 

The boundary conditions on u on £ translate in the following way to boundary 
conditions on u on R n . 

• The Dirichlet condition u = <p on S is equivalent to the Dirichlet condition 
u = <p on R n , where <p := <p o p g L 2 (R n ; C m ). 

• The Dirichlet regularity condition Vsm = <p on S ( Vs denoting the tangential 
gradient on S), is equivalent to V x u = <p on R n , where <p := Po(<p) G 
L 2 (R n ;C nm ). 

• The Neumann condition (V, AVqm) = <p on S (contrary to tradition, v being 
the inward unit normal vector field on S) is equivalent to (eo, AV ttX u) = <p 
on R n , where <p := \ J(p )\(p o p G L 2 (R n ; C m ). 

In this way the Dirichlet /regularity/Neumann problem with coefficients A in the 
Lipschitz domain Q is equivalent to the Dirichlet /regularity /Neumann problem with 
coefficients A in the half space R+ +n , and it is straightforward to extend Theorems 12 .31 
and 12.41 to the Lipschitz domain Q. 



3. Road map to the proofs 

For the reader's convenience, we give in this section an informal explanation of 
the main ideas behind the methods and the proofs of Theorems 12.31 and 12.41 In 
particular, the precise definitions of classes of solutions will be given later. Our 
basic idea for constructing solutions u to the divergence form equation §B) in R^_ +n 
is to consider it as a first order system with the gradient V^m as the unknown 
function. In fact, solving for the t-derivatives in the equation, the divergence form 
equation for u becomes a vector-valued ODE 

d t (V t , x u) + T A (V t , x u) = 0, 

where Ta is an operator only involving the first order derivatives along R n and 
multiplication by entries of A(t,x). The divergence form equation was first studied 
through this ODE in [7]. However, it turns out that if one instead of Vt,x u takes 
the conormal gradient Va« defined by fl3J) as the unknown, then the corresponding 
operator T4 has a simpler structure; the ODE reads 



d t f + DBf = 0, with / := V A u, 



where D 



and B is a second strictly accretive coefficient matrix 



div a 
-V x 

determined by A. This was the key discovery in [9] when coefficients A do not 
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depend on t. This carries over to t-dependent coefficients and as this result is 
central to us, we give full proof of it in Proposition 14.11 

3.1. The Neumann and regularity problems. The first order approach is most 
natural for solving the Neumann and regularity BVPs, since these boundary condi- 
tions are conditions on the conormal gradient /, not on the potential u. Indeed, the 
Neumann BVP means that the normal part (fo)± = d VA u\-^n is given at the boundary 
t — 0, whereas the regularity condition is that the tangential part (/o)n — ^x u \n. n is 
given. Note that both BVPs mean that "one half" of the function / is prescribed. 
This is natural for a first order elliptic equation. 

On the other hand, the set of all traces fo = /|r« of solutions / to the equa- 
tion d t f + DBf = in R+ + ™, with appropriate estimates, is a subspace of L 2 = 
L 2 {H n ; C^ 1+n - )m ) which we denote E\l-L. Here the reader should have the classical 
situation in mind, where (|SJ) is the Cauchy-Riemann equations and E\H is the up- 
per Hardy subspace of L 2 . Just like in this classical situation, it turns out that for 
t-independent coefficients and small perturbations thereof, E\% is a closed proper 
subspace, being "one half" of L 2 , and there is a bounded Hardy type projection E\ 
onto E\H. Moreover, there is a Cauchy type reproducing formula for the solution / 
to the ODE, given f G E\H. Hence there is a one-to-one correspondence between 
solutions / and their traces fo G E\'K. 

Proving these facts for small perturbations of t-independent coefficients is the 
main work in this paper. For t-independent coefficients, this result is in [9]. Before 
explaining the proofs in more detail, assume for a moment the stated properties of 
E\ and E~^H, in order to explain the implications for BVPs. The unique solvability 
of the Neumann BVP means that for each boundary data tp G L 2 (R n ; C m ) there is 
a unique / G E^H such that f ± = ip. Hence well-posedness of the Neumann BVP 
is equivalent to 

E+U^L 2 (R n ;C m ):f^f ± 

being an isomorphism. Similarly well-posedness of the regularity BVP is equivalent 
to 

E\U -+ {g G L 2 (R"; C nm ) ; cni\ x g = 0} : / ^ /„ 

being an isomorphism. Even for t-independent coefficients, these maps are not 
always invertible. For counter examples based on [HE] in this context, see [T2] . 
However, we mentioned three important classes of t-independent coefficients where 
techniques are available to prove invertibility. 

Concerning this approach to solving BVPs, it should be noted that the problem 
is divided into two parts: (i) understanding representation formulas and the trace 
space E\% for solutions to the differential equation, and (ii) understanding the 
relation between the full trace space E\% and the boundary condition (Neumann 
or regularity). Carefully note that (i) only involves the differential equation and not 
the boundary conditions. This is of great importance, since it means that it suffices 
to study the ODE (jHJ), and explains why the trace spaces E\l-i may be well behaved 
even when the BVPs fail to be well-posed. The main harmonic analysis work goes 
into proving that the projections E\ are bounded for all complex t-independent 
coefficients, and small i-dependent perturbations thereof. From this it is deduced 
that the projections E\ depend continuously on A in a certain Carleson sense, and 
as a consequence in Corollary 18.61 that well-posedness of the Neumann and regularity 
problems is stable under small perturbations of the coefficients. 
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We next explain our methods for solving the ODE, i.e. solving problem (i). 
For this we study (jSJ), where we first consider t- independent coefficients B = B , 
and we write E\ = Eq. In this case, we view DB as an unbounded operator in 
L 2 (R n ; C < - 1+n ' )?n ), and at a first glance the solution to (jSJ) with initial/boundary data 
/o seems to be f t = e~ tDB °f . However, the problem is that DB is not a sectorial 
operator, but instead bisectorial, i.e. its spectrum is contained in a double sector 
around the real axis. This indefiniteness means that the operators e~ tDB ° are not 
well defined on L2(R n ; C ( - 1+n ^ m ) for any t ^ 0. Another technical problem is that 
DBq has an infinite dimensional null space. The fact is that there are topological 
splittings 

L 2 = U® U(DB ) = (e+U © EoTij © N(DB ), 

where H = R(DBq) = R(D) is the closure of the range. The splitting of H into 
the spectral subspace Eq % for the sector in the right half plane and the spectral 
subspace Eq"H for the sector in the left half plane is a deep result, and builds on the 
Kato square root problem as discussed in the introduction. This proof also shows 
that DBq has square function estimates, which in particular shows that —DB 
generates a bounded holomorphic semigroup in EqW, and that DB generates a 
bounded holomorphic semigroup in EqH. 

Given any / G Eq"H, differentiation as seen in [9] shows that the generalized 
Cauchy reproducing formula 

/ = Co /0) 

with C^f (t,x) = (e~ tDB °EQ f )(x), yields a solution to (jHJ with trace / . Con- 
versely, given a solution / to (jS|), it is a fact that f t belongs to the range H for any 
conormal gradient / and t > 0. (Note that it follows from (jSJ) that d t f t G %.) We 
apply the projections E^ and suitable exponentials to the equation, giving 

W°E+f a ) = 0, sG(0,t), 
£ > DB °Eof a ) = 0, sG(t,oo). 

Integration with limits lim^oo f s = and lims^o f s = f indicates that the trace f 
belongs to the subspace Eq % and that the generalized Cauchy reproducing formula 
/ = Cq fo holds. This converse and existence of the limits are shown in Corollary 
18.41 For further details concerning the relation between Cq and the classical Cauchy 
integral, see Ex. 1.2] and pi Thm. 2.1]. Note that C + / = for all /„ G E^U. 

For t-dependent perturbations B of a given t-independent coefficient matrix Bq, 
we write the ODE as d t f + DB f t = D£ t f t , with S := B — B . The above argument 
now gives 

■ ds{e _ {t _ s)DBoE + fs) = e -(t-s)DBo E + D£sf ^ s e (Q)t)) 

9s ( e (s-t)DB 0E - f ^ = e (s-t)DBo E - Dg j sj g e (t; ^ 

and integration and subtraction of the equations give the integral equation 

(9) / = C +/i + + S A f, 

for some h + G Eq"H (which will be shown to be h + = EqJq) and where S A is the 
integral operator given by 

i*t /*oo 

(10) S A f t := / e-^ DBo E+DS s f s ds - / e^ DB ° E^ DS.f.ds, t > 0. 

Jo Jt 
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To construct solutions to OH]), it is therefore natural to think of the Picard fixed 
point theorem. For this, we need an appropriate function space of functions in R^ +n 
that contains the free evolution g = C^h + and on which Sa is bounded. By [9], the 
non-tangential maximal function in Definition 12. II of g belongs to L 2 (R n ). Thus the 
space X defined by ||Af*(/)|| 2 < oo is a natural candidate (see Section E]). 

We turn to a closer look at Sa, and that allows us to justify the equivalence 
between ([9]) and the ODE. This operator involves forward and backward maximal 
regularity operators. Usual treatment of maximal regularity is not sufficient for our 
needs. A convenient way to study Sa, in the spirit of functional calculus, is as 
follows. Define, for A G C with Re A > 0, the bounded integral operator -F(A) : 
L 2 (R+; C) -)• L 2 (R + ; C) : u = (u s ) s>0 ^ F(X)u by 

F{X)u t := [ Xe- {t ~ s)x u s ds. 



o 



By letting -F(A) act pointwise in x G R n , it defines a bounded operator on L 2 (R + ; %). 
In this space we also have the operator 



\DB \f:-. 



DB f, f G E+U, 
■DB f, feE H, 



by letting it act pointwise in t G R + . The operator |-D-Bq| is a sectorial operator in 
"H, hence in L 2 (R + ; "H), and on the sector containing its spectrum, A H- F(X) defines 
an operator-valued holomorphic function. Similar to Dunford functional calculus, 
we can apply (F(A))a to \DB \, since \DB \ commutes with each -F(A), and we 
obtain an operator F(\DBq\) on L 2 (R + ; "H). An advantage of this method, referred 
to as operational calculus, is that boundedness of the singular integral / = (f s )s>o l— 
F(\DB \)f defined by for all t > 0, 

F(\DB \)ft ■= f \DB \e-^ DB °\f s ds, 



can be easily derived from the square function estimates for DB , in exactly the same 
way as the boundedness of Ej^, e~ t){DBo \ or more general operators in the functional 
calculus of DB is proved. 

The upshot of all this is a representation 

(11) S A = F(\DB \)E+S + F*(\DB Q \)E^S, 

where Eq are bounded operators such that EqD = (DB )Eq and F*(X)u t := 
J 00 Ae~^ _ ^ A M s c/s. Since all operators on the right hand side in (TTTT) are bounded 
on L 2 (R + ;'H), so is S^. More generally, this boundedness proof goes through on 
weighted spaces L 2 (R + , t a dt; T-L) for \a\ < 1. Details are in Section [HI However, 
for a = — 1, which is the natural scaling for BVPs with L 2 data, boundedness of 
F*(\DB \) fails. Recall nonetheless that the free evolution belongs to X, which con- 
tains L 2 (R + , dt/t; L 2 ) as a subspace. (As compared to the space X, L 2 (R + , dt/t; L 2 ) 
consists of functions that must vanish in a certain sense at t = 0.) Thus we can 
replace L 2 (R + , dt/t; L 2 ) by the larger space X. To ensure boundedness of Sa on X, 
we still need additional control on £ such as finiteness of ||£||c of Definition 12.11 (in 
fact a possibly smaller quantity suffices) and our fundamental estimate is 

(12) \\N.{SjJ)\\ 2 <\\A-Mo\\N.U)h> 
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obtained from the chain 

X A L 2 (R+, dt/t; L 2 ) -% L 2 (R+, dt/t; H) F *^f ol) X. 

See Lemma 15.51 for this modified Carleson embedding theorem in the first arrow 
and Theorem 16.81 for remaining details. This allows us to prove that the trace of 
/ at t = exists in a certain sense and one sees that EqJq = h + . Details of this 
representation of solutions / are in Theorem 18.21 

Now, smallness of this Carleson norm implies smallness of in which 

case (EJ) rewrites as a Cauchy reproducing formula 

/ = (/ - 5 A )~ 1 C + / i + , 

by inverting I — Sa in X. Conversely, for any h + G E^"H, the Picard iteration 
scheme applies to produce a solution f — (I — Sa)~ 1 Cq h + to (Q whose trace / is 
determined by a linear relation E\h + = f . 

This provides us with the Hardy type projections E\ needed to solve the Neumann 
and regularity problems as described above. Details are in Corollary 18.61 Additional 
a "priori square function estimates on solutions, under further regularity assumption 
on the coefficients, can be shown. They are proved in Section 110.21 

3.2. The Dirichlet problem. For the Dirichlet problem with L 2 boundary data, it 
is not obvious that the above first order approach applies. Nevertheless it is possible 
to adapt the arguments and we describe this now. Instead of (jHj) for the conormal 
gradient /, we want to work with the potential u at the L 2 level. The heuristic is 
that u solves the divergence form equation if and only if there is a vector-valued 
potential v such that 

(13) d t v + BDv = 

and u = —v ± . (The minus sign is just for convenience.) On one hand, we have 
seen that u is a divergence form equation if and only if (dt + DB)(V am) = 0. On 
the other hand, applying D to the ODE fl33D gives (dt + DB)(Dv) = 0. Thus 
the claim amounts to rewrite the conormal gradient / = Va« as / = Dv. (Note 
that the latter equation is pointwise in t.) Evaluating the tangential part f n shows 
the relation u = —v ± . Once again the reader should compare this to the classical 
situation of Cauchy-Riemann's equations, in which case u would be a harmonic 
function, and v the analytic function having u and its harmonic conjugate function 
as imaginary and real parts. 

We mention that the v we construct via the conormal gradient again does not quite 
satisfy (1T31) because of the null space of D, but that is not a problem as we only 
need its normal part which does not depend on this null space. Having this heuristic 
in mind we proceed as before in two steps: (i) understanding representation of 
solutions, and (ii) understanding its trace and the relation to solvability. We mostly 
concentrate on (i) as (ii) will follow right away. 

For t-independent coefficients B = B , we note that B D is another bisectorial 
operator, just like DB . It is not injective (the null space is H^) and the spectral 
projections E for the sectors in the right and left half planes split the range BqH 
into two closed subspaces. Similar to the argument for DBq, we have that solutions 
obeying a certain square function estimate to the ODE (1T31) all are of the form 
v = CqV + c for a unique v G EqL 2 and some c G Q( l+n ) m ^ where Cq is defined by 
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(Co~v )(t,x) = (e tB ° D 'EqV )(x). Thus, for t-independent coefficients, we have the 
representation 

u = c - (C+v ) ± , v G £ + L 2 , c G C m , 

for solutions it to (pQ) obeying a square function estimate. This is in Corollary 19.41 
and improves [§] where this was shown for a smaller class of solutions. Note that 
the conormal gradient of u can be calculated as / = DCqVq. 

For t-dependent perturbations B of a given t-independent coefficient matrix Bq, 
suppose we are given a solution u. Since we do not know v yet, we go via / = V A u 
(which in the end will be Dv) to solve ( fl3i) . Since / satisfies (dt + DB)f = 0, as in 
Section 13.11 we look for a functional setting in which we obtain an equation of the 
form ft = e~ tDB °h + + SU/* for some h + in a positive spectral space so as to again 
apply the Picard fixed point theorem to construct solutions. The main difference 
here is in the free evolution term g t = e~ tDB °h + . Indeed, /o, which we should relate 
to h + , may only be defined in a space of Sobolev type with regularity index -1 (and 
in fact, we do not care about fo). 

Square function estimates (see [9]) lead us to the solution. Indeed, we have 
f™\\e- tDBo h + \\ltdt « \\D- l h + \\l provided D~ l h + e L 2 . Here D~ l is defined as 
a closed operator as we assumed h + in a positive spectral space. Hence the natu- 
ral (considering the method) space of conormal gradients for Dirichlet problems is 
y := L 2 (R + ,tdt; L 2 ). Indeed, imposing the free evolution g to be in y allows us to 
obtain g = DCqH + for h + e EqL 2 determined by h + = D~ l h + . (Note that D and 
Bo are interchanged.) This is consistent with the t- independent case. 

The next step is thus to bound Sa in 3^- We use again the operational calculus 
representation (iTTj) . In the space y, it is now the operator F(\DB \) that fails 
to be bounded, but again the additional Carleson control ensures boundedness of 
H^Hy-^y and our second fundamental estimate reads 

(14) \\S A f\W { UfM) A || C ||/|U 2(t(i t;L 2 ). 

Details are in Proposition 17.11 This allows us to justify the formal manipulations 
and to obtain an a priori representation of conormal gradients / = V^w in y of 
solutions u by 

(15) / = DCfh + + S A f, ~h + e E+L 2 . 

We next want to exhibit the announced vector-valued potential v which must 
satisfy / = Dv. Remark that D being non-injective, there is some freedom in the 
choice but we basically want to factor out D in (115j) . This is granted for the free 
evolution and for Sa we obtain starting from (fTUl) that Sa/ = DSa/, where 

/t poo 
e -(ts)B D E + gsfsds _ ^ e (s-t)BoD E - £sfsds 

All this follows from the intertwining b(DBo)D = Db(BoD) of the functional calculi 
of DB and B D. Thus we can set 

v : = Cfh + + S A f, 

and then show that there exists a constant c G C m such that 
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for solutions u to (JT]) obeying an initial square function estimate with / being the 
conormal gradient of u. Again, this is not the only possible choice for v but any other 
choice has identical normal part and this is what we need to recover u. Moreover, 
this choice has good estimates. Details of this representation of solutions u to the 
divergence form equation are in Theorems 19.21 and 19.31 We note that this proves 
that all solutions with gradient satisfying a square function estimates are, up to 
constants, continuous in t with values in L 2 (R n , C m ). This representation also 
enables us to show existence of the trace vq = lim^o^t in L 2 (R n ; CK 1+n ) m ) and h + 
is determined by h + = EqVq. Most importantly, the representation allows us to 
prove non-tangential maximal function estimates in Section 110.11 We remark they 
are a priori estimates. 

Upon smallness of the Carleson control, these representations above rewrite 

/ = (I - S A )- 1 DC+h + ) h + G E+L 2 , 

v = Cfh + + S A f, 

u = c - (v) ± , 

and, conversely, this can serve (via Picard's iteration scheme to obtain the first 
equation) as an ansatz to produce a solution u starting from a given h + G EqL 2 and 
constant c = 0, and Vq G L 2 is determined by a linear relation v = E\h + . With 
this ansatz, well-posedness of the Dirichlet problem is seen to be equivalent to 

E+L 2 L 2 (R"; C m ) : h + h- u = -(v ) ± 

being an isomorphism. This allows us to prove stability of well-posedness of the 
Dirichlet problem under small perturbations in Corollary 19.51 

4. Integration of the differential equation 

Following [U], we construct solutions u to the divergence form system ([T]), by 
replacing u by its gradient g as the unknown function. Consequently (00) for u is 
replaced by (1171) below for g. Proposition 14.11 reformulates this first order system 
(117j) further, by solving for the t-derivatives, as the vector-valued ODE (118j) for the 
conormal gradient 



/ = Vau = [d VA u, Vxuf, where [a, v] 



for a G C m and v G C nm , and d VA u := (AVt, x u)± denotes the (inward!) conormal 
derivative of u. 

According to the decomposition of m-tuples into normal and tangential parts as 
introduced in Section [21 we split the matrix as 

~A ±± (t,x) A ±ll (t,x) 
A ll± (t,x) A u (t,x) 



A(t,x) 



Note that with our assumption that A be strictly accretive on % for a.e. t > 0, the 
matrix A ±± is invertible. 

Proposition 4.1. The pointwise transformation 

A- 1 -A~}A in 



A^A 



An i A. i \ Awn Aw i A i ]A i 
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is a self-inverse bijective transformation of the set of bounded matrices which are 
strictly accretive on TL. 

For a pair of coefficient matrices A = B and B = A, the pointwise map g h-> / = 
[(Ag)j_, <7||]* gives a one-one correspondence, with inverse g = [(Bf) ± ,f l] ] t , between 



, 2 oc (R 



+ ! 



L 2 (R n ; C (1+n)m )) to the equations 



solutions g € L\ 
(17) 

and solutions f G L 2 oc (R + ; "H) to the generalized Cauchy-Riemann equations 



dw ttX {Ag) = 0, 
curing = 



(18) d t f + DBf = 0, 

where the derivatives are taken in Ri +n distributional sense, and D :- 





-V, 



div 2 




This was proved in [HJ Section 3] for t-independent coefficients. The proof goes 
through without changes for t-dependent coefficients, but for completeness we give 
the proof of this important result. Note that R(-D) = "H. 



Proof. We first look at the correspondence A i— y A — B. Fix an arbitrary t > 
and we write A for A(t, •). From the accretivity of A on "H, it follows that the 
component A ±± 

A±± Aj_ 
/ 



.4 
.4 



is pointwise strictly accretive, hence invertible, and therefore so is 
Thus, multiplication by A is an isomorphism on % and, letting 



/ 

_A_L A\\ 

g EH and / = 



i-r-l 



B = AA is bounded if A is so. We calculate, for any fixed 



Ag, 



Re(BfJ)=Re(AAg,Ag) = Re 



A u 
Re 








9± 




A ±± 


Alii 


Al II. 




.011. 


•> 









/ 




^11 




9± 




9± 




A u 




9\\ 




9\\ 



Re(Ag,g). 



This shows that B = A is strictly accretive if A is so. That A = A is straightforward 
to verify, and this shows that A and A are in one-to-one correspondence. 

Next consider a pair of functions g and / in L2 OC (R + ; L2(R n ; C ( - 1+n ^ m )) such that 
/ = Ag. Equations (1T71) for g are equivalent to 



(19) 



d t (Ag) x + div x (A l|x 5( ± + Am,^,, 
%u - = 0, 
curl^^n = 0. 



0. 



The last equation is equivalent to f t G H. Moreover, using that (Ag) ± = f ± , g n = f n 
and g ± = (Bf) ± = A~]_(f ± — A ±n f n ), the first two equations are seen to be equivalent 
to the equation dtf + DBf = 0. This proves the proposition. □ 

Remark 4.2. In terms of the second order equation where g is a gradient and / is 
the corresponding conormal gradient, the identity Re(Bf,f) = Re(Ag, g) rewrites 
Re(£> V au, Vau) = Re(AVf iX M, Vt, x u) for any appropriate u (not necessarily a solu- 
tion). 
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We now want to construct solutions to (fTHl) . Let us first recall the situation when 
B(t,x) = B (x) does not depend on the t-variable. In this case, we view B as 
a multiplication operator in the boundary function space L 2 (R n ;C( 1+n > m ). Define 
closed and open sectors and double sectors in the complex plane by 

S^ + :={\eC ; |argA| <u}U{0}, S„ := S u+ U (S„ + ) , 
S° + := {A e C ; A ^ 0, | argA| < u}, S° u := S° u+ U (S° u+ ), 

and define the angle of accretivity of B to be 

u := sup | arg( J B /, f)\ < vr/2. 

The method for constructing solutions to the elliptic divergence form system, de- 
veloped in j7J[9], uses holomorphic functional calculus of the infinitesimal generator 
DBq appearing in the ODE (118jl . and the following was proved. 

(i) The operator DB is a closed and densely defined w-bisectorial operator, 
i.e. o-(DB ) C S u , where u is the angle of accretivity of B . Moreover, there 
are resolvent bounds ||(A — DBq)~ 1 \\ < 1/dist (A, S u ) when A ^ S^. 

(ii) The function space splits topologically as 

L 2 (R n ; C (1+n)m ) = H H(DB ), 

and the restriction of DBq to % — R(-D) is a closed, densely defined and 
injective operator with dense range in H, with same estimates on spectrum 
and resolvents as in (i). 

(iii) The operator DBq has a bounded holomorphic functional calculus in %, 
i.e. for each bounded holomorphic function 6(A) on a double sector S°, u < 
v < 7r/2, the operator b(DB ) in % is bounded with estimates 

\\b{DB )\\ u ^ u <\\b\\ Loom . 

For background material on sectorial operators (which is straightforward to adapt 
to bi-sectorial operators) and their holomorphic functional calculi, see [2]. The 
construction of the operators b(DB ) is explained in detail in Section 16.14 in the 
more general case of operational calculus. The two most important functions 6(A) 
here are the following. 

• The characteristic functions x + W an d X~(ty f° r t ne right and left half 
planes, which give the generalized Hardy projections := x ± (DBq). 

• The exponential functions e~*' A ', t > 0, which give the operators e~^ DB °^. 
Here |A| := Asgn(A) and sgn(A) := X + W ~ xW- 

A key result that we make use of frequently, is that the boundedness of the projec- 
tions shows that there is a topological splitting 

(20) H = EqH © EqH 

of H = R(-D) = R(DB ) into complementary closed subspaces EqW := R(E^). 

We also recall the definition of the generalized Cauchy extension Cq from Sec- 
tion [31 

Proposition 4.3. The generalized Cauchy extension 
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of /o G EqT-L gives a solution to d t f + DB f = 0, in the strong sense f G 
C 1 (R + ;L 2 ) fl C°(R + ; D(DS )), with L 2 bounds sup t>0 \\f t || 2 ~ H/0II2 -^2 lim- 
its lim^o A = fo and lim^^ / t = 0. 

Now consider more general t-dependent coefficients B(t, x). Fix some t-independent 
coefficients B , strictly accretive on %. (This B should be thought of as the bound- 
ary trace of B, acting in R^ +n independently of t.) To construct solutions to the 
ODE, we rewrite it as 

(21) d t f + DB f = DSf, where £ t := B - B t . 

However, while d t f + DB f = can be interpreted in the strong sense, ( 1211) will 
be understood in the sense of distributions. The following proposition rewrites this 
equation in integral form. It uses operators Eq, defined as 

(22) Eq := EqB^PboH, 

where Pb h denotes the projection onto BqH in the topological splitting L 2 = BqH® 
T-L 1 - and B$ 1 is the inverse of B$ : H. — > BqH. Beware that Bq 1 is not necessarily 
a multiplication operator and is only defined on the subspace B 7i . Note also that 
unlike Eq, Eq are not projections. 

Proposition 4.4. Iff G L 2 oc (R + ; W) satisfies d t f+DBf = in R+ +n distributional 
sense, then 

- [ t r)' + (s)e-( t -^ DB ^E+f s ds= f r) + (s)DB e-^ DB ^E+£ s f s ds, 
Jo Jo 

/oo /*OG 
^{s)e-^ DB «\E^f s ds = J V 4s)DB e-^ DB ^E^£ s f s ds, 

for allt > and smooth bump functions 7)±(s) > 0, where n + is compactly supported 
in (0,t), and rj_ is compactly supported in (t, 00). 

Proof. In Section 13.11 we showed formally how to integrate the differential equation 
and arrived at ©• To make this rigorously, we proceed as follows. By assumption 

(23) + (D<p s ,B f s )y s = J™(D<p s ,£J s )ds, 

for all G C£°(R}+ n ; C( 1+n ) m ). To prove the identity on (0,t), let O G H be any 
boundary function and define <p s := rj + (s)(e~^ s ^ DB ^E^)*(j)o G Cg°(R+; D(D)). To 
show that we can use this as test function, take 77 G C^°(R n ) with 77 = 1 in a 
neighbourhood of x = and J Rn 77 = 1 and write rj e := e~ n r)(x/e) and 

cf>^(x) := V+ (sHx/R) ( Ve * ((e-K*-^lE o ± )*0o)) (x). 

It is straightforward to verify that d s (f) R ' e — > d s <p and D<f> ,e — > Dip in L 2 (supp?7 + x 
R"; C( 1+ ") m ) when R -» 00, e ->■ 0. 
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From 023|) we obtain 

/\-^( s )( e -C'-s)l^°l£+)*0 o - r, + (s)(DB e-^ DB ^E+)*<f> , f s )ds 
Jo 

\ V+ {s)D{e-^ DB ^E^y^B f s )ds 

^ + {s)D{e-^\ DB »\E+y^£ s f s )ds. 



Since B*D(e-^\ DBo \E+)* = (e-^-^^E^DBo)* = {DB e-^ DB °\E^)* , the 
last two terms on the left hand side cancel. Using that Eq D = EqDPb -h — 
Eq (D Bq) Bq 1 Pb h = DBqEq on the right hand side, we have proved that 

h j\' + ( s )^ s)lDB ° lE ofsds^ = (^ j\ + (s)e-^ DB ^E+D£ s f s ds 

Since this holds for all <po, the (0, t) integral formula follows. The proof for the (t, oo) 
integral formula is similar. □ 

Our goal is to take limits to arrive at an integrated equation. Formally, if we 
let r]± approximate the characteristic functions for (0,t) and (t, oo) respectively, we 
obtain in the limit from Proposition 14.41 that 

Eth ~ e'^E+fo = f DB e^ DB ^E+£ s f s ds, 

Jo 

DB e~^ DB ^E^£J s ds, 

if lim^o ft = fo and \im t ^oo ft = in appropriate sense (and yet to be proved). 
Subtraction yields / = Cq fo + SU/, which we wish to solve as 

(24) / = (I - S^C+fo, 

where the integral operator Sa is 



(25) 



l*t POO 

S A ft = / DB e-^ DB ^E+£J s ds - / DB e^ s ^ DB ^E^£J s ds. 
Jo Jt 



and Cq is the generalized Cauchy integral defined via the semigroup e~^ DB °^. 

The equation f[2"4"l) can also be viewed as a generalized Cauchy integral formula, for 
t-dependent coefficients A, and we shall see that, given any fo e L 2 (R n ;C( 1+n ) m ), 
it constructs a solution f t to the elliptic equation. However, for this one and also 
justification of the limiting arguments, one needs a suitable functional setting we 
now introduce. 



5. Natural function spaces 

It is well known that solutions g to ( TT7I) with L 2 boundary data typically satisfy 
certain square function estimates, as well as non-tangential maximal function esti- 
mates. In this section, we study the basic properties of some natural function spaces 
related to BVPs with L 2 boundary data. 
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Definition 5.1. In R+ +n , define the Banach/Hilbert spaces 

X := {f : B}+ n -> C^ m ; #,(/) G L 2 (R n )}, 

POO 

y:={f: Rl +n C< 1+n > m ; / ||/ t ||*tdt < oo}, 

with the obvious norms. Here iV* denotes the modified non-tangential maximal 
function from Definition O By y* = L 2 (R+ +n , dt/t; C^ 1+n ^ m ) we denote the dual 
space of y, relative to L 2 (R+ +n ; & l+n ^ m ). 

In Sections El and M we demonstrate that the maximal function space X is the 
natural space to solve the Neumann and regularity problems in, whereas y is natural 
for the Dirichlet problem. Natural is meant with respect to the method. That the 
spaces y and X are relevant for L 2 BVPs with t-independent coefficients is clear 
from the following theorem. For proofs, we refer to [HI Proposition 2.3] and 
Proposition 2.56]. 

Theorem 5.2. Let fo belong to the spectral subspace E^"H. Then its generalized 
Cauchy extension f = Cq fo as in Proposition \4-3\ has estimates 

\Wh ™ \\f\\* ™ WMW 

We will show in Corollary 18.41 that any distributional solution / 6 X to dtf + 
DB f = is the generalized Cauchy extension of some fo 6 Eq"H. 

Clearly y C L2° C (R+; L 2 ). The following lemma shows that X is locally L 2 inside 
R+ +n as well, and is quite close to y*. 

Lemma 5.3. There are estimates 



1 f 2t ~ f°° ds 

sup - / \\f s \\lds<\\N*(f)\\l< / \\f s \\l-. 
t>o t Jt Jo s 

In particular y* C X . 

Proof. The second inequality follows by integrating the pointwise estimate 

t>0 J JW(t,x) S J J\y-x\<c c lS S 

For the lower bound on ||A^(/)|| 2) it suffices to estimate t -1 J^ ' ||/ s || 2 rfs, uniformly 
for t > 0. To this end, split R" = [J k Qk, where Qk all are disjoint cubes with 
diagonal lengths c\t. Then 

* _1 / / l/M| 2 ^<|Q fc | inf |iV*(/)(*)| 2 < / |iV,(/)(x)| 2 rf.x. 

Jt JQ k x ^ JQ k 

Summation over k gives the stated estimate. □ 

The space y* is a subspace of X of functions with zero trace at the boundary R n , 
in the square L 2 -Dini sense lim^o t~ l J*. 2 * ||/ s |||(is = 0. A fundamental quantity is 
the norm of multiplication operators mapping X into y*. 

Definition 5.4. For functions 8 : R+ +rt -> £(C (1+n)m ), we denote 

:= ||£||#-o>' = sup 
Il/ll*=i 

the norm of pointwise multiplication by 8. 
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The following lemma gives a sufficient Carleson condition for a multiplication 
operator to map into this subspace. 

Lemma 5.5. For functions £ : R^ 4 ""' — > C(C^ 1+n ^ m ) , we have estimates 

ll^lloo < ||£||* < ||£||c> 

where \\£\\c denotes the modified Carleson norm from Definition \2.1\ 

Proof. For the first estimate, fix t and consider only / supported on (t, 2t) in the 
definition of 1 1 ^ 1 1 «v— >-y * - Lemma [5.31 shows that 

sup \\£f\\ y ./\\f\\ x « sMt-^Sfhyit-^fh) = esssup \\£ S \U 

t<s<2t 

where the first two suprema are over all ^ f £ L 2 ((t, 2t) x R n ; C( 1+n ) m ). Taking 
supremum over t shows the estimate ||£||oo ^ 11^- II*- 
For the second estimate, we calculate 

\£f\\y^ II . [jL II dsdy]\£(t,x)f(t,x)\^ 

\£(t,x)f(t,x)\^\ dsdy 




W(s,y) 



±- II \f(t,x)\ 2 dtdx]dsdy<\\£\\ 2 c 

JJW( S ,y) J 



2 

X ■ 



'W(s,y) 

where the final estimate is by Carleson's theorem. □ 

We have not been able to prove that the || • ||* norm is equivalent to the modified 
Carleson norm, that is to prove the appropriate lower bound. It is however easy 
to see that the || ■ ||* norm dominates the standard Carleson norm || • || c . Indeed, 
choosing / as the characteristic function for the Carleson box (0, l(Q)) x Q (times 
a unit vector field) in the estimate < II^IMI/IU, shows that 

||£|| c :=su P -?- // \S{t^)\^ <\\£\\l 

Q M J J(0,l(Q))xQ 1 

Furthermore, it is straightforward to see that the modified Carleson norm is domi- 
nated by the corresponding modified square Dini norm 

f°° dt 

£\\% S, / sup H£J|L — . 



\C ~ / SU P II--II00 . 

q t<S<C()t 

Proof of Lemma \2.2\ We shall prove the lemma assuming only finiteness of the stan- 
dard Carleson norm (hence the lemma also holds for the star norm). Thus we assume 
\\A — A \\ c < oo and do the rest of the proof replacing || • ||c by || ■ || c . 
To prove uniqueness, we use \\A' — A) lie < \\A — A' \\ c + \\A — A || c < oo to obtain 

for all cubes Q C R n , which only is possible if A' = A almost everywhere, by 
t-independence. 
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Similarly for A — A , we have J* (|<2| 1 Jq \A(t,x) — A (x)\ 2 dx)dt/t < oo for any 
cube Q, and it follows that we have essential infimum 



essinf — — / \A(t, x) — An(x)\ 2 dx = 0. 
xt<i(Q) \Q\J fy 



0<t<l(Q) \Q\ Jq 

To prove || v4 II oo < ll^lloo, assume e > and pick Q such that (IQ)" 1 Jq lA^dx) 1 / 2 > 
\\AqWoo - e. Then choose t G (0,l(Q)) such that \Q\~ X L \M^ x ) ~ A (x)\ 2 dx < e 2 . 
We assume that t is a Lebesgue point of t h-> |C}| _1 Jq |A(t, x) — A (x)\ 2 dx and of 
t^lQ]- 1 / Q |A(t,x)| 2 ^. This yields WA^-e < (\Q\~ l $ Q \A Q \ 2 dx) 1 ' 2 < WA^ + e. 
Letting e — >■ proves the claim. 

To prove k > k, assume e > and pick / G H n C °°(R n ; C( 1+n ) m ). Pick Q such 
that supp f C Q. Then 

Re(A /, /) = Re(A/, /) + Re((A - A*)/, /) > 

for all t G (0, l(Q)). Taking the essential supremum over such t gives Re(Ao/, /) > 
kII/IH- Since 7/nC£°(R n ; C^ 1+n ^ m ) is dense in H, taking infimum over /, this proves 
the claim. □ 



6. HOLOMORPHIC OPERATIONAL CALCULUS 

Throughout this section A denotes a closed, densely defined w-sectorial operator 
in an arbitrary Hilbert space "H, i.e. cr(A) C S u +, and we assume resolvent bounds 
|| (A — A) _1 ||^^ < 1/dist (A, S u+ ). For simplicity, we assume throughout that A is 
injective, and therefore has dense range. In our applications A will be |D-B |, and 
% will be the Hilbert space from ©. See Section [3 

The goal in this section is to develop the theory needed to make rigorous the limit- 
ing argument following Proposition 14.41 To this end, we study uniform boundedness 
and convergence of model operators 

(26) St ft := f vt(t, s)Ae-^ A f s ds, 

Jo 

V ;(t,s)Ae-^ A f s ds, 

acting on functions ft(x) = fit, x) in a Hilbert space L 2 (R+, dfi(t); %). For uniform 
boundedness issues, it suffices that the bump functions r/+(t, s) and rj~(t,s) are 
uniformly bounded and compactly supported within {(s, t) ; < s < t} and {(s, t) ; 
< t < s} respectively. For convergence issues and to link to the ODE, they should 
approximate the characteristic functions of the above sets. A convenient choice 
which we shall use systematically is the following. Define rf{t) to be the piecewise 
linear continuous function with support [1, oo), which equals 1 on (2, oo) and is linear 
on (1,2). Then let rj e {t) := r/°(i/e)(l - r]°{2et)) and 

ri?(t,s) :=r l °(±(t-s)/e)r le (t)r le (s). 
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We study the operators from the point of view of operational calculus. This 
means for example that we view = F(A) as obtained from the underlying oper- 
ator A (acting horizontally, i.e. in the variable x) by applying the operator- valued 
function A i— > F(X), where 

(F(X)f) t := [ t V t(t,s)Xe-^ x f s ds, 
Jo 

which depends holomorphically on A in a sector S° + containing the spectrum of A. 
Note that each of these vertically acting, i.e. acting in the t- variable, operators 
F(X) commute with A. 

6.1. Operational calculus in Hilbert space. Consider A as above. Let K. := 
L 2 (R+, d{j,(t); %) for some Borel measure \x. We extend the resolvents (A — A) _1 G 
£("H), A ^ S u+ , to bounded operators on K (and we use the same notation, letting 
((A — A)~ 1 /) t := (A — A) _1 (/() for all / G K and a.e. t > 0). These extensions 
of the resolvents to /C clearly inherit the bounds from %. We may think of them 
as being the resolvents of an cu-sectorial operator A = A^, although this extended 
unbounded operator A^ is not needed below. 
Define the commutant of A to be 

A' := {T G C{K) ; (A - A)~ X T = T(A - A)' 1 for A £ S^}. 

Fix u> < v < 7r/2, and consider classes of operator- valued holomorphic functions 

{holomorphic F : S° + — > A'}, 

{FeH(S: + ;A') ; ||F(A)|| <min(|A|MA|- a ),somea>0}, 
{F G H(S° + ; A') ; sup ||F(A)|| < oo}. 

Through Dunford calculus, we define for F G ^(S° + ; A') the operator 
(28) ^(A):=^ / F(A)(A-A)-MA, 

where 7 is the unbounded contour {re ±l9 ; r > 0}, u < 9 < u, parametrized counter 
clockwise around S^ + . This yields a bounded operator F(A), since the bounds on 
F and the resolvents guarantee that the integral converges absolutely. 

Remark 6.1. Functional calculus of the operator A is a special case of this op- 
erational calculus (128jl . Applying a scalar holomorphic function /(A) to A with 
functional calculus is the same as applying the operator-valued holomorphic func- 
tion F(X) = f(\)I to A with operational calculus. For the functional calculus, we 
write ^(S° + ) and Hqo^S",) for the corresponding classes of scalar symbol functions. 

We also remark that a more general functional and operational calculi for bisecto- 
rial operators like DB are developed entirely similar to those of sectorial operators 
A, replacing the sector by the bisector S u . 

The following three propositions contain all the theory of operational calculus 
that we need. To be self-contained and illustrate their simplicity, we give full proofs, 
although the propositions are proved in exactly the same way as for functional 
calculus, and can be found in [2]. 



H(S: + ;A'): = 
-f^oo(5'° + ; A') := 
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Proposition 6.2. IfF,G G ^(S° + ; A'), then 

F(A)C7(A) = (FG) (A). 

Note that we need not assume that F(X) and commute for any A, fi G S^. 

Proof. We use contours 71 and 72, with angles to < d\ < 62 < 7r/2, so that 72 
encircles 71. Cauchy's theorem now yields 



(W) G (A)=(/^)(/^) 



71 " / \^72 



, 71 A - A \,/ 72 /i — A / ,/ 72 \J 71 11 — A J 11- A 

= [ P^-2mG(\)d\ - = (27n) 2 (FG)(A), 
J 11 A — A 

using the resolvent equation. □ 

Proposition 6.3. Assume that A satisfies square function estimates, i.e. assume 
that 

f°° dt 

/ ||^(tA)w||^— « for all ue V. 

Jo t 

and some fixed ip G ^/(5'° + ). TTien i/iere exists C < 00 swc/i i/ia£ 

||F(A)|| < C7 sup ||F(A)||, / ra«FG^;A'). 

We remark that if square function estimates for A hold with one such ip, then 
they hold for any non-zero ip G ty(S° + ). 

Proof. Note that the square function estimates extend to u G K,, with || • \\jc instead 
of || ■ ||%. We drop /C in || ■ 1 1 /c - Using the resolution of identity J °° ^ 2 {sA)uds / s = cu, 
where < c < 00 is a constant, and the square function estimates, we calculate 



2 dt 



\\F(A)u\\ 2 n I \\ip(tA)F(A)u u 







{ip(tA)F(A)ij}(sA))(il}(sA)u 



ds 







dt 

T 



s 

V(t/s)-) ( I™ V (t/s)M(sA)u» 2 



s° + Jo \Jo s J \Jo s J t 

r°° dt 
<sup||F(A)|| 2 / ||V'(^A)n|| 2 -<sup||F(A)|| 2 ||n| 
s° + Jo s s° 



We have used the estimate 



\mA)F(A)iP(sA)\\ < / ||F(A)|||^(tA)^( S A)A- 1 rfA| < sup \\F(\)\\ri(t/i 
where rj(x) := min{x a , x~ a }(l + | logx|) for some a > 0. 



□ 
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Proposition 6.4. Assume that A satisfies square function estimates as in Propo- 
sition UTc\ Let F n £ (S° + ; A') , n = 1,2, . . satisfy sup n A ||F n (A)|| < oo, and let 
F £ iy oo (S'° + ; A'). Assume that for each fixed v £ /C and A £ we /iave strong 
convergence lim^oo ||F n (A)v — F , (A)v|| = 0. Then the operators F n (A) converge 
strongly to a bounded operator F(A), i.e. 

F n (A)u — >■ F(A)u, for all u £ ]C, as n -> oo. 

Proof. Since sup n ||F n (A)|| < oo by Proposition ^. 31 it suffices to consider u = ip(A)v 
for some fixed ip £ ty(S° + ) \ {0}, since R(^(A)) is dense in K. From (12*51) . we get 

||F n (A)u-F m (A)u|| < / || (F„(A) - F m (A))w|| |^(A)A" 1 rfA|, 



where ||(F„(A) — F m (X))v\\ < ||w|| and |^(A)|/|A| is integrable. The dominated 
convergence theorem applies and proves the proposition. □ 

Propositions 16.21 16.31 and 16.41 show that we have a continuous Banach algebra 
homomorphism 

fU^°+;A')^/;(/C):F^F(A), 

provided that A satisfies square function estimates as in Proposition 16. 3[ This is the 
operational calculus that we need. Note that with F(A) defined in this way for all 
F £ H^S^; A'), Proposition E31 continues to hold for any F n £ if QO (S'° + ; A'). 

6.2. Maximal regularity estimates. Here, we apply the operational calculus from 
Section [6~T1 to prove weighted bounds on the operators Sf 1 from ("""""") and (j27]) . 

Theorem 6.5. The operators are uniformly bounded and converge strongly as 
e — > on the weighted space L 2 (t a dt;TL) if a < 1. The operators S~ are uniformly 
bounded and converge strongly as e — > on the weighted space L 2 (t a dt; Ti) if a > — 1. 

Note that the case a = is the usual maximal regularity result in L 2 (dt; Ti). The 
methods here provide a proof of it. 

To establish boundedness of the integral operators F(X), we rely on the follow- 
ing version of Schur's lemma. The proof is straightforward using Cauchy-Schwarz' 
inequality. 

Lemma 6.6. Consider the integral operator f t i— > J °° k(t, s)f s ds, with C-valued 
kernel k(t,s). If the kernel has the bounds 



/'OC poo 

— \k(t,s)\s Pl ds = Ci < oo, sup^— / \k(t,s)\t p2 dt = C 2 < oo, 
2 a Jo s s l:>1 a J 

for some fli, f3 2 £ R, then the integral operator is bounded on L 2 (t a dt; Ti) with norm 
at most \/C\C 2 . 

The second result that we need shows that when the integral operators F(X) define 
a holomorphic function in ^/(S 1 ^; C(JC)), then the resulting operator F(A) can be 
represented as an integral operator with operator-valued kernel. 

Lemma 6.7. Consider a family of integral operators F(X)f t = J °° k\(t, s)f s ds such 
that the C-valued kernels have the bounds 

sup^y o Mt,s)\s^ds<rj(X), su P^y o Mt,s)\t^dt< V (X). 
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// sup AgS o 77(A) < 00, if A !->■ k\(t,s) is holomorphic in S° + for a.e. (t,s), and 
if ff K \d\k\(t, s)\dtds is locally bounded in A for each compact set K, then F £ 
H x (Sl + \C(L 2 (rdt\U))). 

If furthermore 77(A) < min(|A| a , |A|~ a ) /or A £ 5'°, and some a > 0, t/^en F £ 
$(,S° + ;/:(L 2 (t a rft; %))), and 

POO 

F(A)f t = / k A (t,s)f s ds, for all f £ L 2 {t a dt]U) and a.e. t, 
Jo 

where the operator-valued kernels k\(t, s) are defined through ( TJ%|) for a.e. (t, s). 

Proof. Schur's lemma |6"1)1 provides the bounds on F(X). To show that the operator- 
valued function F is holomorphic, by local boundedness it suffices to show that the 
scalar function 

A^ JJ(h t ,k x (t,s)f s )dsdt 

is holomorphic, for all bounded and compactly supported /, h. The hypothesis on 
d\k\(t, s) guarantees this. 

To prove the representation formula for F(A), it suffices to show that for each 
/ £ L2(t a dt;'H), v £ H, and a.e. t, changing order of integration is possible in 



(v } k x (t, s)(X - A) _1 /s)dsdA. 
Thus, by Fubini, one needs to show 



J J \k x (t,s)\\\f s \\ds 



\d\\ 

< 00, tor a.e. t. 



|A| 

The bounds on k\(t, s) in the hypothesis guarantee this. □ 

Proof of Theorem \6.5[ Since in L 2 (t a rft;'H) and S~ in L 2 (t _Q ctt; %), with A re- 
placed by A*, are adjoint operators, it suffices to consider S+. Let 

F e (X)f t := f\t(t,s)Xe-^ x f s ds. 
Jo 

Uniform boundedness of the integral operators F e (X) follows from Lemma [6.61 with 
0i = —a, 02 — 0, using the estimate J Q e x x~ a dx < e y y~ a , which holds if and only if 
a £ (—00, 1). Indeed, since A £ S° + with v < tt/2, we have Ai := Re A ~ |A| and 

/ \Xe~ x(t ~ s) \s~ a ds « f X 1 e~ Xl{t - s) s- a ds = X«e~ Xlt [ 1 e x x~ a dx <t~ a , 
Jo Jo Jo 

Similarly, \Xe- x ^\dt < e AlS / A ~ e~ x dx = 1. 

Again using Lemma [6.61 we note for fixed e > the crude estimate ||F e (A)|| < 
|A|e- eReA , and with Lemma EH we verify that F e £ ^(S° + ; £(L 2 (t a dt;H))), and 

F £ (A)f t = [ r)t(t,s)Ae-^ A f s ds = St ft, for a.e. t. 
Jo 

To prove strong convergence, by Proposition 16.41 it suffices to show strong con- 
vergence of F e (X) as e — > 0, for fixed A £ S° + . By uniform boundedness of F e (X), 
it suffices to show that F t (X)f converges in L 2 (t a c?t; "H) as e — > for each / in 
the dense set U<5>o -^((^ t a dt; H). This will follow from norm convergence of 
F e (X) in C(L2((5,5- 1 ),t a dt;n),L 2 (t a dt;n)) for each fixed 5 > 0. To see this, we 
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use Lemma [6.61 with j3i = —a and 02 = 0. As above C\ is uniformly bounded. One 
verifies decay to as e — > of 

poo ps+2t 

sup / X ie ~^ s)Xl dt and sup / X^^^dt. 

This shows that C2 — > as e — > 0, which proves the strong convergence and the 
theorem. □ 

6.3. Endpoint cases. The operators S~ are not uniformly bounded on L 2 (£°<i£; ~H) 
when a < —1, and therefore no limit operator S~ exists in these spaces. Indeed, if 
r)(t) is a smooth approximation to the Dirac delta at t — 1 and / G "H, then S~(r)f) t 
is independent of e for e < £/2, with non-zero value ~ Ae _A / G H for £ »j 0. Thus 
su Pe>o Jo°° \\^e(vf)t\\nt a dt = 00 if a < — 1. By duality cannot be uniformly 
bounded when a > 1. 

In this section we study the endpoint cases a = ±1. It is convenient here (and as 
we apply these abstract results in the subsequent paper [6]) to introduce the abstract 
spaces Y := L 2 (tdt;U) and Y* := L 2 (dt/t;%). Note that they differ from y, y* 
by the target space H being here an arbitrary Hilbert space. To obtain a uniform 
boundedness result for S~ , assume there exists an auxiliary Banach space X with 
continuous embeddings 

(29) Y* CX C L 2 oc (dt;U)i 

i.e. fa\\ft\\udt < \\f\\x ^ Jo°° ll/tllw^A no ^ 101 eacn nxe ^ < a < 6 < c», and 
such that the map u H- (e~* A M)( >0 is bounded % — > X, i.e. 

(30) l|e~ <A «|U ^ WAu, for all ueE 



Theorem 6.8. Consider the model operators and S~ from WbXWl) and Y , Y* 
and X as above. 

The operators are uniformly bounded on Y* and converge strongly to a limit 
operator S + G C(Y*, Y*) as e ->■ 0. 

The operators S~ are uniformly bounded Y* — > X , and there is a limit operator 
S~ G C(Y*,X) such that lim e ^ \\S~f — S~ fW^a^-H) — for any fixed < a < b < 
00 and f eY*. 

For the proof, we shall need the first part of the following lemma. The second 
part will be required in Propositions 17.11 and 17.21 below. 

Lemma 6.9. The operators 

POO 

/ r] e (s)Ae~ sA f s ds :Y* -*U 
Jo 

are bounded, uniformly in e, and converge strongly as e —> 0. Let U s : % — > Ti 
be bounded operators such that H —> Y*;h !->■ (U*e~ sA * h) s>0 is bounded. Then the 
operators 



00 



r] e (s)e- sA U s f s ds:Y 
are bounded, uniformly in e, and converge strongly as e — > 0. 
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Proof. For the first operator, square function estimates for A* give 

ds 



For the second operator 



H 



sup 

IW|2 = 1 



(sA*e~ sA hJ s ) Ve (s)- 



< 



TkfWv < ll/lh 



i$ sup 

H l|fc||2 = l 



(U:e- sA *h,f s )r) e (s)ds 
< sup ||*7; e - A */i|| K .|MI 



Y < \\Vef\\ 



Y < 



5 • 



where in the second last estimate the hypothesis is used. (Note that the %-bound 
on U s is not used quantitatively.) 

To see the strong convergence, replace rj e by rj t — r\ e i and use the dominated 
convergence theorem. □ 

Proof of Theorem \6.S\ The result for is contained in Theorem 16. 5 \ so it suffices 
to consider S~ . Write 



(31) S;f t 



V ;(t,s)Ae-^ A f s ds 

»t+2e 



7/7 (t, s)A(e- {s - t)A - e- {s+t)A )f s ds 
( Ve (t) Ve (s)- V ;(t,s))Ae-^ A f s ds 

POO 

+ Ve(t)e~ tA / r] e (s)Ae- sA f s ds = 1-11 + III. 
Jo 



We show that it is only the last term which is singular in the sense that it is 
not uniformly bounded on Y*. Consider the term / and the symbol F^(X)u t = 
r)~(t, s)k\(t, s)u s ds, where k\(t, s) := Ae~^~ t)A (l —e~ 2tx ). Boundedness of Fj(X) 
on Y*, uniformly in A and e follows from Lemma l6~6l and the estimates \k\(t, s)\sds < 
t and f Q s \ k\(t, s)\dt < 1. For example 



\k x (t,s)\sds < min(l,tAi)e tAl J Aie _sAl sds = tmin(l,Ui)(l + l/(tAi)) < t, 

with Ai := Re A. On the other hand, for fixed e > 0, it is straightforward to 
verify with Lemma [6.61 that \\F t :(\)\\y*-*y* ^ |A|e _eReA , and with Lemma 16.71 that 
Fl e V(S°;£(Y*)) and 



rj~(t, s)A(e-^ A - e^ s+t)A )f s ds, 



for a.e. t. 



To prove strong convergence, as in the proof of Theorem l6.5[ by uniform boundedness 
it suffices to show norm convergence of -F/(A) in jC(L 2 ((5, S" 1 ), t~ 1 dt] H), Y*) for each 
fixed 5 > 0. This follows from Lemma 16.61 where one verifies decay to as e — > 



of sup^^-i) J Q e \k\(t, s)\dt and sup^^-i) f*_ 2e \k\(t, s)\dt, and hence C 2 -+ 0, for 
fixed A G S° + . Together with the uniform bound sup t t~ l f°° \k\(t, s)\sds < oo, this 
proves the strong convergence for the term I. 
Consider next the term II and the symbol 

rt+2e 

F e n (X)u t = / (Ve(t)Ve(s) - v7(t, s))Xe-^ x u s ds. 
Jo 
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Boundedness of -F/ 7 (A) on Y*, uniformly in A and e follows from Lemma EISl and the 
estimates \\e~( s+t ^ x \sds < t and f°?„\Xe~^ s+t ^ x \dt < I. On the other hand, for 

JO 1 1 ~ Js/3 I I ~ ' 

fixed e > 0, we verify with Lemma [6761 that ||F/ / (A)||y.^y. < |A|e _eRcA , and with 
Lemma E2] that F e H e ^(S° + ;£(F*)) and 

F e H (A)f t = / (Ve(t) V e(s) - V7(t, s))Ae-^ A u s ds, for a.e. t. 
Jo 

With the same technique as for the term /, the strong convergence of the term // 
follows from the decay to as e — > of sup^^-i) J^_ 2e \\e~( s+t ^ x \dt. 

It remains to estimate the principal term III. Since the variables t and s separate, 
we can factor this term through the boundary space T-L as a composition Y* — > H — > 
X, where Lemma I6~9l and the assumed bounds e~* A : "H — > X prove boundedness, 
uniform in e, as well as strong convergence as maps Y* — > H — > -^2(0? b; H). This 
completes the proof. □ 

7. Estimates of the integral operators Sa and Sa 

Let us come back to our concrete situation. Consider the operator DB from 
SectionH We set A = \DB \ := DB sgn(DB ) on "H = R(-D), and see that A satisfies 
the assumptions of Section^ It is a closed, densely defined, injective operator with 
cr(A) C and ||(A — A)" 1 ^^ < 1/dist (A, S u+ ) (this follows from the resolvent 
bounds on DB ). We apply the abstract theory from Section l6~3l to this A and spaces 
Y* := y* n 4 oc (R+;%), X := X n 4° C (R+;H) and Y := y n 4° C (R+;H). Note 
that the continuous embeddings (1291) follow from Lemma 15.31 and the boundedness 
hypothesis (1301) on H — > X : h h-» (e~*' DS °'/i) f> o follows from Theorem 15.21 (and the 
analogous result for the lower half space Ri +n , i.e. fo e EqT-L giving a solution of 
d t f + DB f = for t < 0.). We shall use the operational calculus of A to rigorously 
define and estimate the operator Sa in (|2*5|) . 

The strategy for the Dirichlet problem described in Section l3~2l leads us to consider 
the functional calculus of BqD and the integral operator Sa from ( fl6i) . If Bq were 
invertible on all L2, then DBq and -Bo-D would be similar operators, but this is not 
the case in general. Still, whenever B Q is strictly accretive on H, it is true that B D 
is an w-bisectorial operator with resolvent bounds. Furthermore, the L 2 space splits 
as 

L 2 = B n © u 1 - 

(cf. (122jl ) and BqD restricts to an injective closed operator with dense range in BqH. 
This operator has square function estimates, and therefore bounded functional and 
operational calculus in BqH, as in Section 16.11 For proofs and further details, see 
[8]. We set A := \B D\ and Eq := x ± (-Bo-D)- We extend an operator b(B D) in the 
functional calculus to an operator on all L2 by letting b(BoD) = on T-L 1 - = N(BqD). 
With this notation Eq(BqH) = and we shall prefer the latter to ease the 

notation. 

A important relation between the functional calculus of DB and B D is 
(32) B b{DB ) = b(B D)B , 

where we also extend operators b(DB ) to all L 2 , letting b(DB )\ N ( DBa ) := 0. The 
equation ( 132]) clearly holds for resolvents b(z) = (A — z)" 1 . The general case follows 
from Dunford integration (|28|) and taking strong limits as in Proposition ^. 41 (adapted 
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to bisectorial operators). Note that (1321) in particular shows that for appropriate b 
and u 

b(DB )Du = Db(B D)u. 

A final observation is that with A = |-D-B | and A = \B D\, then A* = \DBq\ and 
A* = \BqD\. So A and A* are of the same type, and the same holds for A and A*. 
The boundedness result for the operator Sa is as follows. 

Proposition 7.1. Assume that £ : R+ +n — > £(C^ 1+n ^ m ) satisfies \\£\\* < oo, and 

define operators 

pt POO 

S A ff-= / r l t^s)Ke-^ K E+£ s f s ds+ / V ;(t, s)A e -^ A E^£ s f s ds. 
Jo Jt 

Then \\S A \\x^x < ||£||* and \\S A \\y^y < ||£||*, uniformly for e > 0. In the space X 
there is a limit operator Sa = S* G C(X; X) such that 



iimH^/ - S A f\\ L2 ^ b . L2 ) = 0, for any f G X,0 < a < b < oo. 

In the space y, there is a limit operator Sa = S A G C(y; y) such that 
Km\\S e A f-S%f\\y = 0, for any fey. 

Proof. The result on X is a direct consequence of Theorem I6.8[ since ||£ f\\y* < 
\\£\\*\\f\\x and S A = S?E+£ + S~E^£. Note that R(E$) C % C L 2 . 

Consider now the space y. The second term S~Eq£ is bounded on y, uniformly in 
e, and converges strongly on y. This follows from Theorem 16.51 and the boundedness 
H^^lly^y < Halloo ^ ||£||* < oo. The term S^Eq£ we split as 

f vt(t, s)Ae-^ A E+£ s f s ds = f v +(t, s)A( e -^ A - e~^ A )E+ £ s f s ds 
Jo Jo 

poo 

- / ( Ve (t)Ve(s) - T}+(t, s))Ae-^ A E+£J s ds 

Jt-2e 

poo 

+ Ve (t)Ae~ tA / Ve (s)e- sA E+£ s f s ds. 
Jo 

The result for the first two terms follows from the proof of Theorem 16.81 by duality, 
only using the boundedness of £ on y. For the last term, as the variables t and s 
split, it suffices to show uniform boundedness and convergence of 

L 2 -> y : h r] e (t)Ae- tA h 



and 



POO 

y^L 2 :f t ^ / Ve (s)e- sA E+£ s f s ds 
Jo 



separately. For the first operator, this follows directly from the square function esti- 
mates for A. To handle the second, it suffices to estimate B J °° r) e (s)e~ sA EQ £ s f s ds = 
Jo°° Ve( s ) e ~ sA EQ£ s fsds, since B is strictly accretive on % D R(c~ sA Eq). To this end, 
we apply Lemma RHU with U s := Eq£ s P h , where P% is orthogonal projection onto 
H, and A replaced by A. The hypothesis there on boundedness of 

H^Y* U* s e- sA *h = P n £* s e~ slDB » l x + {DB*)h, 
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follows from the maximal estimate in Theorem 15.21 (with B replaced by Bq), the 
assumed boundedness of £* : X — > y* and L 2 boundedness of x + (DBq) and 
This completes the proof. □ 

By inspection of the proofs above, the limit operator Sa, both for / G X and 
/ G y, is seen to be 

S A ft = lim (j^ Ae- (t - s)A E+£ s f s ds + jf Ae-^ A E^£ s f s ds 

with convergence in L 2 (a,b; L 2 ) for any < a < 6 < oo. This holds since we may 
equally well choose to work with the characteristic function rj°(t) = X(i,oo)(t) instead 
of the piecewise linear function 77° defined below ( 12611271) . The only places we need 
the continuity of 77° are in Theorems 18.21 and 19.21 below. 

We now turn to defining rigorously the second integral operator needed for the 
Dirichlet problem. Write Cb(X, V) for the space of bounded and continuous func- 
tions on X with values in V. 

Proposition 7.2. The operators 

S\ff.= / vt(tiS)e-^ A E+£J s ds- / rj;(t,s)e-^ A E^£ s f s ds 
Jo Jt 

are bounded y — > Cb(R+;L 2 ), with sup t>0 H^/tl^ < ll^ll*ll/l|y> uniformly for e > 
0, and there is a limit operator Sa G £(y,Cb(H + ;L 2 )) such that \im t ^ \\S A f t — 
SU/tlh = locally uniformly for t G (0,oo) ; for any f G y. The limit operator 
satisfies Sa/ = DSa} in distributional sense, where Sa = Sa is the operator 

from Proposition \7.1\ and has limits 

lim \\S A f t -h~\\ 2 = 0= lim ||SU/t||2, 

t— >0 t—>oo 

where h~ := — / °° e~ sA E ~£ s f s ds G Eq L 2 , for any f G y. 

Note that S A ft = when t £ (e, e _1 ), so convergence S A ft — > SAft is not uniform 
up to t = 0. By inspection of the proof below, the limit operator is seen to be 

rt poo 

(33) S A ft = / e-^ A E+£J s ds - / e~^ A E^£ s f s ds, 

Jo Jt 

where the integrals are weakly convergent in L 2 for all / G y and t > 0. 

Proof. The estimates for Sa are more straightforward than those for Sa since there 
is no singularity at s — t. For the (0, t)-integral, split it as 

f nt(t, s)e^ K (I - e~ 2s ~ A )E+£ s f s ds + e~ tA [ n+{t, s)e~ sA E+£ s f s ds. 
Jo Jo 

For the first term, we write e-^ A (I-e- 2sA ) = T ^- s ((t-s)Ae- ( - t -^ A )((I-e- 2sA )/(sA)) 

to obtain the estimate \\e^^~ s ^ A (I — e~ 2sA )|| < s/t. From this uniform boundedness 
and convergence, locally uniformly in t, as e — > follows by Cauchy-Schwarz in- 
equality. For the second term we use uniform boundedness of e~ tA and duality to 
estimate it by 

i-t 



sup 

W\2=l 



[ (£:e- sK \E+yh,f s ) V Ut,s)ds 
Jo 



< \\S* 
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using Lemma 16.91 as in the proof of Proposition 17. 11 Moreover, the L 2 difference 
between the integral at e and e' is bounded by f* \\f s \\ 2 \vt(t, s) — r)£(t, s)\ 2 sds — > 
as e, e' — > for fixed t, which proves the convergence. 

The proof for the (t, oo)-integral in S A is similar, splitting it as 

/oo _ /'oo 

V ;(t, s)e-^' A (I - e~ 2tA )E^£ s f s ds + e~ tA J ^(t, s)e~ s ' A E^ £J s ds, 

and using the estimate ||e _ ^~*- )A (/ — e _2tA )|| <t/s for the first term and Lemma f6T9l 
for the second. 

Since clearly S A f G C;,(R + ; L 2 ), its locally uniform limit S A f also belongs to 
C&(R + ; L 2 ). To find the limits of S A f t at and oo, since S A : y — > l7&(R + ;L 2 ) 
is bounded it suffices to consider / G y such that f t — for t ^ (a, b), with 
< a < b < oo fixed but arbitrary. In this case, 

S A ft = I e~^ A E+£ s f s ds - I e~^ A E^£J s ds 

J a<s<min(t,b) >/max(f ,a)<s<6 

satisfies E^SAft = for t < a and Eq SAft = when t > b, from which the two 
limits lim t ^o Eq SAft = = lim^oo EqSa^ follow. For the remaining two limits 
lim^so E$S A f t and lim^ Eq S A f t , we use that 

lim / \\e^ t - s)A E+£J s \\ 2 ds = Q = \im [ ||(e" (s -' )A - e~ sA )E^£J s \\ 2 ds 

t— >oo I t— >0 I 
J a J a 



by dominated convergence. 

To verify the identity S A = DS A , note that / °° (<p t , S A f t )dt = J °° (D(f> t , S A f t )dt 
for all / G y and G C£°(R^ +re ; C^ +n ^ m ). This uses the relations BqE$ = 

B q E±Bq 1 P BoH = EoPb q h = Ef by (JS2D- Let e ^ and use S\ and S A con- 
vergence. This completes the proof. □ 

8. The Neumann and regularity problems 

Throughout this section, A denotes t-dependent coefficients satisfying ^ and ©, 
and Aq G Loo(R n ; £(C^ 1+n ^ m )) denotes t-independent coefficients which are strictly 
accretive on T-L. We let B := A and Bq := Aq be the transformed strictly accretive 
coefficients from Proposition 14.11 and define £ := Bq — B . 

For the Neumann and regularity problems, one seeks estimates of the gradient 
9 = V^w rather than the potential u. With a slight abuse of notation, we say 
below that g solves the divergence form equation when u does so. 

Definition 8.1. By an X-solution to the divergence form equation, with coefficients 
A, we mean a function g G L l 2 oc (~R + ; L 2 (R n ; C (1+n)m )), with estimate ||iV*(#)|| 2 < oo, 
which satisfies 

dw t , x Ag = 0, 
curing = 0, 

in R^_ +n distributional sense. 

Note that the boundary behaviour of g is not specified in this definition; we show 
existence of a limit in appropriate sense (see also Section [TT|) . This will allow us 
to formulate in what sense the boundary data is prescribed. Our representation 
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and the boundary behavior of ^-solutions bears on the following result. Recall that 
A = \DB \. 

Theorem 8.2. Assume that < oo. Let f G X. Then f G L^ C (R + ;H) and 

solves d t f + DBf = in R+ +n distributional sense if and only if f satisfies the 
equation 

f t = e ~ tA h + + S A f t , for some h + G E+H. 
In this case, f has limits 

p2t i-2t 

(34) hint- 1 / - fogds = = lim r 1 / \\f s \\ 2 2 ds, 

where f := h + + h~ and h~ := J °° Ae~ sA Eg £ s f s ds G E^Ti, with estimates 

max(||/i + || 2l ||/|-|| a )w||/o|| 2 <||/|U. 
If furthermore ||£||* is sufficiently small, then there are estimates 

\\h-h<\\h + h^\\foh^\\f\\x. 

Proof, (i) Assume that / G X n L 1 2 oc (R + ;'H) and satisfies the differential equation. 
To show that ft = e~ tA h + +S A ft, we choose rjf for rf 1 in Proposition H3 and subtract 
the equations to obtain 

/t POO 
(d sV t)(t,s)e-^ A E+f s ds + J (d sV ;)(t,s)e-^ A E^f s ds 

rt poo 

= / V ?(t,s)Ae-^ A E+£ s f s ds+ / V ;(t,s)Ae-^ A E^£J s ds. 
Jo Jt 

Note that DB = ±\DB \ = ±A on EqH. WefixO<a<6<oo and consider the 
equation in L 2 (a, b; %). By Proposition 17.11 the right hand side converges to S A f in 
L 2 (a, b; H). When t G (a, b) and e is small, the left hand side equals 

(36) e" 1 J* e~ sA (E+f t _ s + E~ f t+s )ds 

- e~ l ^ e~^ A E+f s ds - 2e f e~^ A E^f s ds. 

Je J^)- 1 

To prove that the first term converges to / in L 2 (a,b;'H), adding and subtracting 
the term e~ x J^ e e~ sA ftds = e~ eA (eA) _1 (J — e~ eA )ft shows that the square of the 
L 2 (a, b; H) norm of the difference is bounded by 

/b / j _ g -eA\ 2 rb r2e 

[I - e~ €A eA j ft dt + J^ 1 J e \\ft-E+f t _ s -Euf t+s \\ 2 2 dsdt^O 

as e — > 0, using Proposition 16.41 for the functional calculus, dominated convergence, 
and the identity f t = Eq f t + Eq f t as / is H- valued. 

Next consider the last term in (1361) . For any <f> G L 2 (a, b; H), we have 



fief 6 e-^ A E^f s ds,(j>t)dt 

Ja \ J(2e)- 1 J 

= e f (f s , f\e-^ A ' - e" sA *)(E ")>^t + e^ A *(E ")* f <p t dt] ds. 

J(2e)~ 1 V Ja Ja J 
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From the sup — L2 estimate in Lemma l5~3l for /, the estimate ||e — ( - s— *-* A * — e" sA * || <t/s 
and the strong limit lim^oo e~ sA * (Eq)* = 0, it follows that the last term in (I36j) 
converges weakly to 0. Hence the middle term must converge weakly in ^2(0, b; L2) 
as well, and we may replace e _ ^ _s ^ A by e~ tA since ||e _ (* _s ) A — e * A 1 1 < s/t. We get 
that 

j* (e-^e- 1 jf * E+f s ds), <f^j dt = (V 1 jT E+f s ds, J" e~ tA * fad^J 

converges for all G L2(a,b; Li\ Since e _1 J 2e Eq f s ds are uniformly bounded in 
% by Lemma and since functions e~ tA * <J) t dt are dense in BqH pa H* (for 
example J^ e e~ tA e~ l <j)dt — > P BoH <j)), it follows that e _1 J^ e Eq f s ds converges weakly 
to a function h + G Eq"H, and that the weak limit of the middle term in (IHtjj) is 
e~ tA h + . In total, this proves that f t — e~ tA h + = SAft- 

(ii) Conversely, assume that / G X and f t = e~ tA h + + SAft- First, / is H-valued 
since e~ tA h + G Eq"H and SAft £ H for almost every t. To verify that / satisfies the 
differential equation, note that (dt + DBo)e~ tA h + = 0. It suffices to show that for 
<p G C™(R\ +n ; C( 1+ ") m ) there is convergence 



J (-d t <j> t + B*D<p t , ft)dt -> J (D<f> s , £J s )ds, 



e ->■ 0, 



where / t e := S^/t. For the term S^E^Ef, Fubini's theorem and integration by parts 
give 

^00 i*t 

/ / V t(t,s)((-dt + A*)^Ae-^ A E^£ s f s )dsdt 
Jo Jo 

/OO / POO \ 
s)d t {e-^ A "^(t> t )dt, E+Esf.j ds 

/OO / POO \ 
U {d t r 1 t)(t,s)e^ t - s ^A* ( l )t dt 1 E+£ s f s j ds 

POO POO 

-> / (A> s ,£ + £J,)cfe = / (D(f> s ,E+£ s f s )d S . 
Jo Jo 

Adding the corresponding limit for the term S~Eo~£f gives the stated result. Note 
that E+ + Eq = P BaH and DP BoH = D. 

(hi) To show the limits, note that Eq f - e~ tA h + = S + E^£f G y*, and by 
inspection of the proof of Theorem l6.8l we see that Eq f — e~ tA J °° Kc~ sA Eq £ s f s ds G 
y*. From this, the limits for / follow. 

(iv) It remains to prove the estimates. Note that ( |20|) and Lemma 15.31 show that 



/2t 
\\fs\\lds< 



2t 

2 

X- 



Proposition 17. II shows that HSUIU^A" < 1/2 if \\£\\* is sufficiently small. In this case 
/ — Sa is an isomorphism on X with ||(J — Sa)~ x \\x-*x < 2. Using this together 
with Theorem 15.21 we get estimates \\f\\x — ~~ SA)~ 1 e~ tA h + \\ x pa \\h + \\2- This 
proves the stated estimates and completes the proof. □ 
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Theorem 8.3. Assume that ||£||* < oo. Then g is an X -solution to the divergence 
form equation with coefficients A if and only if the corresponding conormal gradient 
f = [(Ag) ± , </||]* G X satisfies the equation 

f t = e ~ tA h + + S A ft, for some h + G E+U. 
In this case, g has limits 

p2t p2t 

(37) limr 1 / \\g s - g \\ 2 2 ds = = lim r l / \\g s \\ 2 2 ds, 

where g := [(-B /o)x, (/o)nP and \\g \\ 2 < \\g\\x holds. If furthermore is suffi- 

ciently small, then there are estimates 

\\h + h ~ \\goh ~ \\g\\x- 

Note that these limits are stronger than L 2 convergence of Cesaro means t~ l J 2t g s ds 
(see Section [TT]). 



Proof. The equivalence follows right away from Proposition 14.11 and the first part of 
Theorem 18. 21 Note that ||<7t||2 ~ \\ft\\2 and \\g\\x ~ \\f\\x- The limit for g at t = oo 
is immediate from that of /, so is the limit of the tangential part g n of g. The limit 
of g x at t — follows from that of Bf. To see this, write 

B t ft - B f = B (f t - fo) - S t f t . 

Since £ t ft G y*, we have lim^o^" 1 L * ||^s/s||2^ s — 0. The limit of i?o(/t — /o) at 
t = follows from the limit of /. The rest of the proof is immediate. □ 

We note the following immediate corollary to Theorem 18.31 

Corollary 8.4. Assume that coefficients A = A are t-independent. Then g is an 
X-solution to the divergence form equation if and only if the associated conormal 
gradient f is a generalized Cauchy extension C^h + of some h + G E^T-L, i.e. 

f t = e ~ tA h + for all t > 0. 

In this case, h + = lim^o ft in L 2 sense. In particular, the class of X -solutions 
in Definition \8.1\ coincides with the class of solutions in 0, Definition 2.1(i-ii)] for 
t-independent coefficients. 

That the solutions in [U] are of this form was shown in the proof of [HI Theorem 
2.3]. Note that the operator Ta\h used in [2j is similar to our operator DB \ n , as 
in Definition 3.1]. 

Remark 8.5. We may ask whether y* could be used as a solution space for gradients 
of solutions. The answer is no because we have seen that functions in 3^* vanish in 
some sense at the boundary so that free evolutions e~ tA fo G y* if and only if fo = 0. 
A second question is then how far gradients of solutions are from being in 3^*- 
Inspection of (iii) in the proof of Theorem 18.21 reveals that ft — e~' A /o G y*, i.e. the 
free evolution e~ tA fo G X is the only term responsible for / (hence g) to belong to 
X and not to y*. 
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We are now ready study BVPs. We recall that for the class of solutions used in [5], 
with t-independent coefficients A Q , well-posedness of the Neumann and regularity 
problems was shown to be equivalent to the maps 



being isomorphisms respectively. From Corollary I8.4[ it is equivalent to well-posedness 
in the class of ^-solutions. 

We now turn to t-dependent perturbations of the coefficients. 

Corollary 8.6. Assume that the Neumann problem for Aq is well-posed. Then there 
exists e > such that for any t-dependent coefficient matrix A with \\£\\* < e, the 
Neumann problem is well-posed for A in the following sense. 

Given any function if G L2(R n ; C m ), there is a unique X-solution g to the di- 
vergence form equation with coefficients A, whose trace g satisfies (A g ) ± = if. 
Moreover, this solution has estimates 



The same holds true when the Neumann problem is replaced by the regularity problem 
and the boundary condition (A g )j_ = (f is replaced by (go)\\ = f £ L 2 (R n ;C nm ) 
such that curLj^ = 0. 

Proof. Throughout the proof, we assume that ||£||* is small enough, so that I — Sa 
is invertible on X by Proposition 17. f I To solve the Neumann problem, we compute 
/ by making the ansatz 



using Theorems l8.2l and l8.3l We see that / satisfies the Neumann boundary condition 
(/o)± = f ii and only if h + solves the equation T^/i" 1 " = if, where Ta '■ EqU — > 
L 2 (R n ; C m ) is the operator 



with / given by fl38|) . Since TA h + = (h + ) ± , a computation using Lemma [6791 and the 
boundedness of (J — Sa)" 1 on X implies \\Ta — ^a \\l 2 ^L2 ^ ll^-ll*- By assumption 
Ta is an invertible operator, and thus Ta remains an isomorphism whenever ||£||* 
is sufficiently small. Thus, in this case we can, given ip, calculate h + = T^ 1 ^ with 
II h + 1| 2 ~ || f || 2 and find a unique solution g to the Neumann problem, with estimates 

\\g\\x ~ 1Mb ~ \\h + h » 1Mb- 

For the regularity problem, we proceed as for the Neumann problem, but instead 
solve for h + in the equation (h + + J °° ke~ sA EQ 8 s f s ds j = if. □ 



E+U -> L 2 (R n ; C m ) : h + ^ (h + ) 



E+U L 2 (R n ; C nm ) ; crnlj = 0} : h + h- (h + ) 



N *(g)h ~ Ibolh « \\<ph- 
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9. The Dirichlet problem 

Throughout this section, A denotes t-dependent coefficients satisfying and (jHj), 
and Aq G Loo(R n ; C(C^ 1+n ^ m )) denotes t-independent coefficients which are strictly 
accretive on %. We let B := A and Bq := Aq be the transformed strictly accretive 
coefficients from Proposition 14.11 and define S := Bq — B. 

Definition 9.1. By a y-solution to the divergence form equation, with coefficients 
A, we mean a function u G W 2 ' loc (R^ +n , C m ), with estimate / °° \\g t \\ltdt < oo of its 
gradient g := V^w which satisfies div^Ag = in R+ +n distributional sense. 



We will prove in Theorem I9.3[ under ||£||* < oo, that any ^-solution belongs to 
C(R + ; L 2 ), modulo constants. Note also that we do not assume any limits of u at 
t = or t = oo, but will prove such below. This will allow us to formulate in what 
sense the boundary values are prescribed. When discussing C(R + ;L 2 ) limits and 
bounds of 3^-solutions u, we shall identify the function u with a C(R + ; L 2 ) function 
modulo constants. 

Our representation of 3^-solutions bears on the following result. Recall that A = 
\B D\. 

Theorem 9.2. Assume that \\S\\* < oo. Let fey. Then f G L l £ c (R + ;H) and 
solves d t f + DBf = in R + +n distributional sense if and only if f satisfies the 
equation 

f t = De~ tA h + + S A ft, for some h + G E+L 2 . 
In this case, let v t := e~ tA h + + SAft- Then f = Dv, v satisfies the equation 

(39) d t v + BDv = -PSf 

in R+ +n distributional sense, where P = I — Eq — Eq = I — Pb h is the projection 
onto 'H ± along BqT-L, and v t has L 2 limits 

(40) lim \\v t — v Q \\ 2 = 0— lim ||t> f || 2 , 

t— >0 t— >oo 

where Vq := h + + h~ and h~ := — J Q °° e~ sA E ~S s f s ds G Eq L 2 , with estimates 
max(||ft, + || 2 , \\h~\\ 2 ) ps ||v || 2 < sup \\v t \\ 2 < \\f\\ y . 

If furthermore ||£||* is sufficiently small, then there are estimates 

\\h~h < \\h + \\ 2 « sup ||ut|| 2 « \\f\\y. 

Proof, (i) We assume that / G 3^nL 2 oc (R + ; %) and satisfies the differential equation. 
As in the proof of Theorem X6.2\ we aim to take limits e — > in equation flHol) . By 
Proposition 17.11 the right hand side converges in y to S A f. Fix < a < b < oo. 
For t G (a, b) and small e, the left hand side equals 

(41) e" 1 j* e- s \E + f t „ s + Eq f t+s )ds 

- e" 1 r e-^ A E+f s ds - 2e f e~^ A EQ f s ds. 
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As in the proof of Theorem 18.2} the first term converges to / in L 2 (a,b; L 2 ). The 
L 2 -norm of the last term is bounded by ej^e)- 1 ll/slk^ ^ e (/(2e)- 1 ll/s|| 2 s(is) 1//2 , and 
hence converges to 0, uniformly for t G (a, b). 

We conclude that ft := e" 1 J e e~^ A E+f s ds converges in L 2 (a, b; L 2 ) as e — > 0. 
In fact, since sup t>0 ||e~* A ||L 2 -s>L2 < oo we have 

1/2 

^—J \\e- {t °- t)A Ul-ff)hdt<( I \\f t e -f7\\tdt 

when to > b. Hence, since (a, b) is arbitrary, fl converges in L 2 , locally uniformly in 
t. Call the limit / and note that it coincides with / — S A f G y for a.e. t > 0. Fix 

the definition of ft . The estimate 



to > and note that f t+to = lim e _>o e tA ft — e tA ft and that in fact f to G Eq"H by 



/■oo 

sup / ||e- tA / t oll^<ll/lS<ll/lS 

t >0 Jo 



follows. Consider the restriction A + of A to EqH, which is a closed and injective 
operator with dense domain and range. We claim that f to G D(A; 1 ). To see this, 
by duality it suffices to show that 

S ll^lla. for all 0G D((A; 1 )*). 

As in the proof of Proposition E31 we use an identity f Q °°(tA + e~ tA +) 2 f ta dt/t = 4 _1 / 4o 
to estimate 



l((A+)*0,/*o)l 



Jo 1 



£ \\9\\2\\J \\y. 



Hence the claim. As D(A + 1 ) = R(A + ) C R(-D), this shows that f ta = Dhf o , where 
h~i G BqEqI-L = EqL 2 has bounds \\hf \\ 2 < \\f\\y, uniformly in to- From the identity 

ft+to = e- tA f t0 = e- tA D~h+,™ get 

j {foj t+t0 )dt= De~ tA *<p t dt,h^ , 

for any G L 2 (a,b; L 2 ). Here the left hand side converges as t — > 0, and the 
functions De~ tA * <j) t dt are dense in H. (For example f^ e De~ tA e~ 1 (j)dt —> D<p.) 
Since ||^||2 is uniformly bounded, it follows that hf -+ h + weakly in EqL 2 as 
to -> 0. Letting t -+ in f t+to = e~ tA Dht = De~ tA ht o , we obtain f t - S A ft = ft = 
De- fA h + for a.e. t > 0. 

(ii) Conversely, assume that / G y and /t = De~ tA h + + S^/*. The right hand 
side is H- valued, so / G L2° C (R+;'H) as well. As in Theorem 18.21 we verify that / 
satisfies the differential equation, and we omit the details. 

(iii) Introduce v t = e~ tA h + + SAft, so that / = Dv and the stated limits on v 
follow from Propositions 16.41 and 17. 2[ To prove ( 1391) . compute 

d t (e- tA h + + S A f t ) = -B De- tK ~h + - B DS A f t 

»2e 



e- 1 J e- sA (E+£ t _J t _ s + E^S t+s f t+s )ds, 
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for e << t. This uses the result for the operator B D analogous to Proposition 14.31 
The claim follows by letting e — > 0, using Proposition 17. II for convergence of DS A = 
S A , and a calculation as in part (i) of the proof of Theorem 18.21 for the last term. 

To prove the estimates, note that the square function estimates for B$D and the 
accretivity of Bq on "H show that 

fh+\\ 2 » \\B De-H + \\ y « \\De- A h+\\ y < \\f\\y + \\S A f\\ y < \\f\\ y . 

From Proposition 17.21 we also obtain the estimates max(||/i + || 2 , ||/?~|| 2 ) ~ \\v0W2 < 
su Pt>o II II 2 i$ ||^ + ||2+ \\f\\y ^ \\f\\y, where we have used the topological splitting 
BqH = EqL 2 © EqL 2 in the first equivalence. 

(iv) Finally, Proposition 17.11 shows that || jS^Hy— *-y < 1/2 if \\£\\* is sufficiently 
small. In this case I — S A is an isomorphism on y, giving the estimate 

\\f\\y<\\De-H+\\ y . 

As \\De~ tA h + \\ y ps \\h + \\ 2 , this proves the stated estimates and completes the proof. 

□ 

We can now prove a rigidity theorem for 3^-solutions. 

Theorem 9.3. Let u be a y -solution to the divergence form equation, with coeffi- 
cients A. Assume that ||£||* < 00. Then there is a constant c £ C m such that 

u = c — v ± 

almost everywhere, where v £ C(R + ; L 2 (R n , C ( - 1+n ^ m )) is the vector-valued potential 
from Theorem \9.m\ obtained from the conormal gradient f of u. Identifying the func- 
tions u and c — v ± , there are L 2 -limits lim^o \\ut — W0II2 = 0, u := c — (v ) ± , and 
lim^oo \\ut — c|| 2 = 0, and there are bounds 

||«0 ~ C|| 2 < SUp \\u t - c|| 2 < ||V t>a .u||y. 

t>0 

If furthermore \\£\\* is sufficiently small, then with h + as in Theorem \9.S\ 

\\h + \\2 ~ ||Vt, x li||y. 

Proof. Let / = [(AV t , x u)±, V^u]* £ y H L 2 oc (R + ; U) as in Proposition O and then 
let h + and v be as in Theorem 19.21 For the equality u = c — v± , it is enough to 
show that Vi.^M = — Vf jX w x in R^ +n distributional sense. It is clear that V ' x u = 
/11 = — V x f ± . Using (|3"9l . we have (9 t t> x + (BDv) ± = 0, because normal parts of 
functions in H 1 - are zero. Since (BDv) ± = (Bf) ± = (V t;X u)± = d t u, we conclude 
that d t u = —d t v ± . 

The stated limits and bounds now follow from Theorem l9.2l and \\f\\y ~ ||V t)X «||y. 

□ 

The constant c in Theorem 19.31 can also be calculated as the limit 

c = lim (u, T d (p), r d <f)(t, x) := <p{t - d, x), 

d— >oo 

for any £ C^°(R 1+n ; C m ) with J <pdtdx = 1. In particular this limit does not 
depend on 0. So if this limit is zero, we obtain a solution that vanishes at 00 in 
averaged sense. (This is thus equivalent to vanishing at 00 in L 2 sense as defined 
in Section [2j) This is akin to the classical pointwise limit at infinity required to 
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eliminate constants for representing harmonic functions in the upper half-space. 
Here the averages replace the pointwise values. 

We also note the following corollary to Theorems 19.31 and 19.21 

Corollary 9.4. Assume that coefficients A = Aq are t-independent. Then u is a 
y -solution to the divergence form equation if and only if it is the normal part of a 
generalized Cauchy extension C$h + of some h + G EqL 2 , modulo constants, i.e. 

u t = (e~ tA h + )± + c for all t > and some c G C m . 

In particular, the class of y -solutions in Definition \9.1\ that vanish at oo in L 2 
sense coincides with the class of solutions in (9j Definition 2.1 (hi)] for t-independent 
coefficients. 

That the solutions considered in [U] are of this form follows from P, Lemma 4.2] 
and the proof of [HI Theorem 2.3]. Note that the operator T A \n used in [9] is similar 
to our operator B D\b h, as m 0, Definition 3.1]. 

We are now ready to study BVPs. We recall that for the class of solutions used in 
[9], with t-independent coefficients A , well-posedness of the Dirichlet problem was 
shown to be equivalent to the map 

E+L 2 -)• L 2 (R"; C m ) : h + h+ (h + ) ± 

being an isomorphism. From Corollary 19.41 it is equivalent to well-posedness in the 
class of 3^-solutions. Remark that well-posedness implies that the map u y u t is 
a Co-semigroup on L 2 (R n ; C m ) and this corollary also shows that the results in [S] 
concerning the domain of this semi-group obtained for solutions in the sense of [§] 
apply to 3^-solutions. 

We now turn to t-dependent perturbations of the coefficients. 

Corollary 9.5. Assume that the Dirichlet problem for A is well-posed. 

Then there exists e > such that for any t-dependent coefficient matrix A with 
||£||* < e ; the Dirichlet problem is well-posed for A in the following sense. 

Given any function ip G L 2 (R"; C m ), there is a unique y -solution u to the diver- 
gence form equation with coefficients A, with boundary trace Uq = (p. Moreover, this 
solution has estimates 

\\^t,xU\\y w sup ||ut|| 2 w \\ip\\ 2 . 
t>0 

Proof. Throughout the proof, we assume that ||£||* is small enough, so that I — Sa 
is invertible on y by Proposition 17.11 To solve the Dirichlet problem, we make the 
ansatz 

(42) u= ( K (I + S A (I-S A r 1 D)Cfh + ^ 

for some h + G E^L 2 . Theorems 19.21 and 19.31 show that u is a ^-solution to the 
divergence form equation with coefficients A and that all 3^-solutions with L 2 trace 
are of this form. Moreover, the Dirichlet boundary condition u = <p is satisfied if 
and only if h + solves the equation T A h + = <p, where T A '■ EqL 2 — > L 2 (R"; C m ) is 
the operator 

f A : h + ^ (h + - jT e~ sK E^£ s f s ds 
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where /:=(/ — S A)~ 1 DCQh + . Since TA h + = (h + ) ± , Lemma HHH and the bound- 
edness of (I — Sa)' 1 on ^ imply that — r^ || i2 ^x, 2 < \\S ||*. By assumption 
is an invertible operator, and thus Ta remains an isomorphism whenever ||£||* 
is sufficiently small. Thus, in this case we can, given </?, calculate h + = T^ 1 ^ with 
II h + 1| 2 ~ II V 9 II 2 and find a unique solution u to the Dirichlet problem. From Theo- 
rem 19.31 we get estimates 

1Mb = IKIh < sup ||wj 2 < ||V t ^M||y ~ ||/i + ||2 « IM| 2 . 

This proves the theorem. □ 

Remark 9.6. The tangential part V\\ of the vector-valued potential 

v = (I + S A (I-S A )- 1 D)Cfh + 

can be viewed as a set of generalized conjugate functions to the Dirichlet solution 
u. Our proof of Theorem 19.21 above eliminates the need of the technical condition 
at 00 on these conjugate functions which was required in [HI Definition 3.1]. 

10. Further estimates 

In Section [HJ we constructed solutions, with estimates on the modified non- 
tangential maximal function, to the Neumann and regularity problems with L 2 
boundary data, and in Section [9] we constructed solutions, with estimates on the 
square function, to the Dirichlet problem with L 2 boundary data. In this section, we 
prove two theorems which give modified non-tangential maximal function estimates 
for the Dirichlet problem, and, upon some further regularity on the coefficients, 
square function estimates for the Neumann/regularity problems. 

10.1. Maximal function estimates for 3^-solutions. 

Theorem 10.1. Let A : R+ +n -> £(C (1+n)m ) with \\A\\oo < 00 and strictly accretive 
on %, and assume that there exists t-independent coefficients Aq with \\A — Aq\\c < 
00. Then any y -solution u to the divergence form equation with coefficients A, with 
boundary trace Uq G L 2 (R n , C m ), has modified non-tangential maximal estimates 

K||2 < ||iV*(u)|| 2 < ||V t>as u||y. 

The core of the proof reduces to the following estimate of the operator Sa- 
Lemma 10.2. For any fixed p € [1, 2), the operator Sa has estimates 

K((S A f) x )\\ 2 <\\£\\ c \\f\\y. 

Here N%(f)(x) := sup <>0 t~ ( - 1+n ^ p \\f\\ Lp ^ w ^ ttX )) is an L p modified non-tangential max- 
imal function. 

Proof of Theorem \10.1\ modulo Lemma \10.2[ As in Theorems 19.21 and 19.31 any y~ 
solution u with L 2 trace (hence, vanishing at 00 in L2 sense) can be written 

u t = (e- tK h + + S A f t ) ± , h + e E+L 2 , fey. 

From Poincare's inequality \\u - u w ^ x) \\ L2{ w^ x )) % t\\V 8jy u\\ L ^w{t,x)), where u W (t,x) 
denotes the average, we obtain the estimate ||A^(tt)|| 2 < || N}{u) || 2 + || Vt^itHy. 
Theorem 15.2} Lemma 110.21 and Theorem 19.21 now apply to give the estimate 

ll^(«)|| 3 <||fc + || a +||/||y« lively. 
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To see the first estimate, write h + = Boh + with h + £ EqH, and apply Theorem 15.21 
to get \\e~ tA B h + \\ x = \\B e- tA h + \\ x < \\h + \\ 2 « \\h + \\ 2 . The lower estimate follows 
from Lemma 15.31 since 

\\ N *( u )\\l ~ limr 1 / ||u s |||ds = ||tt ||l. 
t^-o j t 

□ 

Proof of Lemma MO . 6 A Before we start, we remark that p y \\N^((S A f) ± )\\ 2 is in- 
creasing, so it suffices to consider p close to 2. We shall fix the value of p eventually 
in (iii) below, when we see how close to 2 it need to be. Next it suffices to prove the 
inequality for t i— > f t compactly supported in R + . Indeed, combining Lemma 15.31 
and Proposition 17.21 for all e > and / £ y we have (since p < 2) 



y- 



< / < lne sup ll^/Hl < lne 11 

Thus, if fs := X^.s- 1 )^)/ f° r / £ 3^ we have for fixed e > 

W(X( e , e -i)(*)(&/) Jlh < liminf \\W(xi e , e -i)(t)(S A f6W\\2. 

Now our assumption gives 

\\NZ( X{ e,e-i)(t)(S A f 5 ) ± )\\ 2 < \\S\\c\\f 5 \\y < \\S\\c\\f\\y, 
uniformly in e, so for all / £ y and e > we obtain 

\\N!(x(^MSAf)dh<\\S\\cU\\y- 

It remains to let e — > and apply the monotone convergence theorem, 
(i) We now fix t (->• f t compactly supported in R + and write 

S A ft = / e-^E+SJJs - / e^ K E^£J s ds =: I - II. 
Jo Jt 

Most of the time we use the pointwise inequality N$ < N*. It is only for one term, 
estimated in (iii) below, that we require p < 2. 
Split the integral I as 

/ = / e~^ K (I - e- 2sK )E+£ s f s ds + e~ tK [ e~ sK E+£ s f s ds = h + I 2 . 
Jo Jo 

As in the proof of Proposition ^. 2^ the kernel of I± has bounds s/t, giving the estimate 

(43) ||Ar. W |g<f ||,<<f (/'^) (/>,,,||^) f 



poo 

< / \\£ s fs\\lsds < \\swU\fWy. 

Jo 



Similarly we split II = Ih+II 2 by writing e"^^ = e-^'^il -e~ 2tA ) + e^ tA e~ sA , 
and a Schur estimate similar to (1431) give the bound for Hi. Next we write 

Ih = e~ tK / e- sK E^£ s f s ds - e~ tA / e~ sk E^£ s f s ds =: II, ~ Ih- 
Jo Jo 
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By Theorem 15.21 the term II 3 has bound 



N* [ Bqc Bq PbqH I e~ sA E^£ s f s ds 



sup 

\h\\2=l 



< 



e~ sA E^£ s f s ds 



(S* s e- sA (Ev)*h,f s )ds 



2 

< \\£\ 



y- 



(ii) It remains to consider I 2 + Ih = (Eq + E )e tA / * e sA £ s f s ds. Note that 
(Eq + Eq) = Pb h- Since we only consider the normal component of I 2 + II4 and 



(P BoH h) ± = h x for any h, it remains to estimate e tA f Q e sA £ s f s ds. To make use of 
off-diagonal estimates (see Lemma [I0.3p . we need to replace e~ tA by the resolvents 
(I + itB D)~ 1 . To this end, define ipt(z) '■= e - *^' — (1 + itz)~ l and split the integral 

t*t _ /*oo i*oo 

e~ tA / e~ sA £ s f s ds = ip t (B D) / e~ s ' A £ s f s ds - / t/; t (B Q D)e- aK S 8 f 8 ds 
Jo Jo Jt 

ft _ /-t 

-l/„-sA t\c c j„ 1 /r 1 ,,-j-d n\-l 



(/ + itB Q D)-\e- sA - /)5 s / s cZs + (/ + zt^oD)^ 1 / S a f a d8. 
'0 Jo 

For the first term, square function estimates show that ip t (B D) : L 2 — > y* C X is 
continuous, and a duality argument like for II 3 gives the bound. For the second and 
third terms, we note the operator estimates 



\\MB D)e 



-sA\ 



t e -t\BoD\ _ (J + itBoD)" 1 



tB Q D 



sB D)e- slBoDl 



< 



t/s, 



and 



\\(I + itB Dr\e- sA -I)\\< 



tBnD e - S \ B ° D \ - I 



1 1 + itB D sBqD 



< 



s/t. 



Schur estimates similar to ( )43i) give the iV* bounds, 
(iii) It remains to prove the estimate 



N*( (I + itBoD)- 1 / £J s ds 



< \\£\ 



cwj \\y- 



To show this, fix a Whitney box W(t ,x ), take h G L g (W(t , x ); C (1+n)m ), and 
let /i = outside W(t ,Xo). Here 1/p + 1/q = 1, p < 2 and q > 2. To bound the 
L p (W(to, #o)) norm, we do the duality argument 



1 rcoto / pt \ 

/ (z+^r 1 / £ s f s d S ,h t )dt 

to Jc-ho V Jo J 



co*o / 1 /"Coto 

- / (/ - itDBlY l h t dt J 

y *0 Jmax(co 1 t ,s) 

< / / l^y)!!/^)!^)^, 



where 



H(y) 



1 /-coto 



Kl-itDB^htiy^dt. 
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To handle the tails of (/ — UDBq) 1 h t , we split the space into annular regions 
R n = Uk=o A k, where A) := B(x ;t ) and A k := (2 k A ) \ {2 k ^A Q ) for k > 1. 
Define f k (s,y) := X(o,c t )(s)XA k (y)f(s,y) and H k (y) := XA k {y)H(y). Then Whitney 
averaging as in the proof of Lemma 15.51 gives 

COt °\£(s,y)\\f(s,y)\H(y)dsdy<f2 [[ \8{s,y)\s\f k {s,y)\H k {y)^- 



t 



(t,x) 

1 f f , „ l0 \ 1 ^ ( 1 /" /" i " |2 ^ 



f^ J JnX+^w^x) \t 1+n J J W (t,x) J V l+n JJw(t,x) 

oo „ 

y^ll^llc / A ( -tttt^II sf k \\ L2 (w(t,x)) -7z+x\\ H k \\L 2 (w(t,x))) {z)dz 



oo 

< 

k=0 

oo 

< 

k=0 



^qj+^W^kW L 2 {W(t,x))j {z)dz 

< ll^llcM^llMR^IIMd^i 2 ) 1 / 2 !!^^). 



k=0 



Here A denotes the area function Ag(z) := {JJ\ y _ z \ <ca |<7( s > y)\ 2 s ^ 1+n ^dsdyY^ 2 and 
N*g(z) := sup| J/ _ 2 | <cs \g(s, x)\ is the non-tangential maximal function, where c G 
(0, oo) is some fixed constant, and M is the Hardy-Littlewood maximal function. 
On the fourth line we used the tent space estimate by Coifman, Meyer and Stein in 
[I7J Theorem 1(a)]. Since M : L q / 2 — > L q / 2 is bounded, we have 



\\M{\H k \ 2 f' 2 \\ Lq{m < \\H\\ Lq{Ak) < - / \\(I-itDB*y%\\ LqiAk) dt 

■y rcoto 

<2 m — \\ht\\L q (B(x ;c t ))dt < 2" m to q \\h\\ Lq (y/(t , xo )). 



The third estimate uses Lemma 110.31 below, and thus is where we choose p < 2 
sufficiently close to 2 so that 2 < q < 2 + 5. We obtain the maximal function 
estimate 



TVf ((I + ttBoD)- 1 [ SJ s ds) (x ) < \\£\\csnpT2- km t n /q - n \\A(sf k )\\ Lp{Il 

V Jo / * 0>0 fclo 

oo , x 1/p 

< ||^|| c ^2- fc ( m -^)sup -—- / \A{sj)\Hx] 

u-n *o>0 \{* to) JB(x :(2 k +cco)t Q ) J 



k=0 *o>0 t J JB(x ;(2 k +cc )t ) 

< \\S\\ c M(A(sfy) 1/p (x ), 
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where c is the constant from the definition of A and m > n/p. Since M : L 2 / p —> L 2 / p 
is bounded, this yields 



<||5|| c ||M(^( a /)P) 1 /i'|| 2 <||£|| c || 1 A( a /)|| a 
This completes the proof of the maximal function estimate. 



\S\ 



c 



y- 



□ 



The following lemma, which we used above, is contained in [TJ Lemma 2.57]. 
However, we give a more direct proof here, since the algebraic setup in [7j was quite 
different. 



Lemma 10.3. Let B be t -independent coefficients, strictly accretive onH = R(D). 
Then for each positive integer m, there is C m < oo and 5 > such that 



'l + itDB )- l f\\ Lq{E) < 



(1 + dist (E,F)/ty 

for allt > and sets E,F C R n such that supp / C F , and all q such that \q—2\ < 5. 
Here dist (E, F) := inf {\x — y\ ; x G E,y G F}. 

Proof. For q = 2, these off-diagonal estimates can be proved as in [HI Proposi- 
tion 5.1], using estimates on commutators with bump functions (and replacing 
the operator B D there by DB ). By interpolation, it suffices to estimate ||(1 + 
2tD-Bo)~ 1 /||L 9 (R n )^L,(R n ); uniformly for t and q in a neighbourhood of 2. To this 
end, assume that (/ + itDB )f = f. As in Proposition 14. H but replacing d t by 
(it) -1 , this equation is equivalent to 

{Aog) ± + itdxv x (Aog)n = (A g)±, 
9w ~ itV x g ± = g», 



where A ,g,g are related to B ,f,f, respectively, as in Proposition 14.11 Using the 
second equation to eliminate g n in the first, shows that g x satisfies the divergence 
form equation 

* \g x =\\ ztdivJ \ A± A ±9± 



Lg ± := [l itdivj A 



By the stability result of Sneiberg jH] it follows that the divergence form operator 
L is an isomorphism L : VK 1 (R ri ) — > W~ 1 (R n ) for \q — 2| < 5, giving us the desired 
estimate 



1 91 



< 



\9i 



t\\V x gj 



\gw\\q ^ \\g\U 



□ 



10.2. Square function estimates for Absolutions under t-regularity for the 
coefficients. Looking closely at the equation div t , x Ag = 0, it seems unlikely that Ab- 
solutions g would in general satisfy the square function estimate J °° ||9^t|||tcZt < oo, 
i.e. d t gt G y, when A is t-dependent. More precisely, it writes d t (A ±± g ± + A ±ll g n ) + 
div x (A n± g ± + A|m<7||) = 0, and as d t and multiplication by A do not commute the 
quantity d t gt does not arise. We show in the next result that d t gt G 3^ can be obtained 
upon a further t-regularity assumption on A. This also improves the regularity of g% 
itself. We do not claim this assumption is sharp nor necessary (in particular, it could 
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well be that this regularity on the components A ±± ,A ±n suffices). This regularity 
assumption is akin to the one in [33], going back to [TH], towards property of the 
L-harmonic measure with respect to surface measure on bounded Lipschitz domains 
for real elliptic operators. See also pU |25] for results with smallness assumptions 
of the derivatives of A. The difference wih these works is that we are imposing our 
coefficients A to be perturbation of "good" t-independent coefficients. So our next 
result neither contains or is contained in the above cited works. Besides, it is again 
an a priori estimate on solutions, so it is valid independently of solvability issues. 

Theorem 10.4. Let A : R+ +n ->■ £(C (1+n)m ) with WA^ < oo and strictly accretive 
on H, and assume that there exists t-independent coefficients A with \\A — A \\* 
sufficiently small. If A satisfies the t-regularity condition 

\\td t A\\* < oo, 

then any X -solution g to the divergence form equation, with coefficients A, with 
boundary trace go has regularity dtgt G L 2 oc (R+j -^2) with estimates 

\\9t9t\\y ^ h\\x- 

We also have estimates sup f>0 ||<7t||2 ~ IMU? and t i— >■ g t G L 2 is continuous with 
limits lim^o \\9t — fi'olh = = lim^oo 1 1 ^ 1 1 2 - The converse estimate \\g\\x % \\dtg\\y 
holds for all X-solutions g, provided \\tdtA\\* is sufficiently small. 

If max(\\tdiA\\*, \\td t A\\*) < 00 for some i — 1, . . . , n, then dig t G L 2 oc (R + ; L 2 ) for 
any X -solution g to the divergence form equation with coefficients A, with estimates 
\\di9t\\y ^ IIpIIat- The estimate \\g\\x ^ ll^xfl'lly holds for all X-solutions g , provided 
||tVt )a; A||* is sufficiently small. 

We do not know whether the smallness assumptions are needed for the converse 
estimates to hold. We also remark that the same conclusions hold for the conormal 
gradient /, as is clear from the proof below. 

Lemma 10.5. If h G X has distribution derivative dth G y, then dt(SAh) G y with 
estimates 

\\dt(S A h)\\ y < (\\S\U + \\td t S\U)\\h\\ x + \\S Hooll^lly- 

Proof of Theorem \10.J\ modulo Lemma \10.5\ (i) As in the proof of Corollary EEJ any 
A'-solution can be written g = [(Bf) ± , /J*, where 

(/ - S A )f = e~ tA h + , for some h + G E+U. 

Introduce the auxiliary Banach space Z :— {h G X ; d t h E y} C X, with norm 
II^IU := \\h\\x + a ll^^lly- By Proposition 17.11 and Lemma [10 .51 we have estimates 
||5U^IU < C||^IU and ||9 t (i5' J 4/i)||y < Z?||/i||^ + C||£^/i||;y, where we assume C < 1, 
and we choose the parameter a > small enough so that 

\\S A \\ z ^z<C + aD<l. 

Hence I — S A is invertible on both X and Z. Since e~ tA h + G Z by Theorem 15.21 
we conclude that / G Z with estimates ||9t/||;y < ll/IU ~ ll e ~ <A ^ + IU ~ ll^lh- F° r 
the gradient g, this gives the bound \\d t g\\y < \\td t B\\J\f\\ x + (H-BIU + < 

\\h + h*\\f\\ X a\\9\\ X . 

(ii) To prove the sup — L 2 estimate and trace result for g t , write J Q °° sr](s)d s g s ds = 

/ °°(?7(s) + sr)'(s))g s ds, for some r] G C£°(R + ). Take the limit as 77 approaches the 
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characteristic function for (0, t) to get 



If* If 1 

gt — - g s ds H — / d s g s sds, a.e. t > 0. 



The last term has bound (L \\dggg\\ 2 sdsy/ 2 , whereas the first term satisfies 



1 

t 



g s ds - 5f 



< 

2 



00 ( i r 2l ~ H \ 



as t — > 0. Hence the trace claims follow from the square function estimates ||9 t ^||y < 
oo. Moreover, the estimate sup t>0 \\gtW2 ^ \\g\\x + ll^fi'lly ^ \\g\\x follows. The 
converse estimate follows from Theorem 18.31 

An integration by parts, similar to above, shows that 



1 f 2t 



2#2t = 9t + ~ J g s ds + - J d s g s sds, a.e. t > 0. 

Taking limsup^^ of both sides, shows 2 lim sup^^ ||gt||2 = umsu Pt->.oo II 5*11 2- Since 
|| g t \\ 2 is bounded, we conclude that lim^oo \\gtW2 = 0. 

(iii) To show \\g\\ x < \\d t g\\y, consider / satisfying e~ tA h + = f t - S A f t . Theo- 
rem 15.21 and Lemma 110.51 give 

\\h+\\ 2 « H^e-^+Hy < \\d t f\\y + + ||tM|.)ll/ll* + \\£\U\d t f\\y, 

where by Theorem 18.21 we have \\f\\x ~ ||/?- + ||2 as is assumed small enough. If 
in addition is sufficiently small, then we obtain \\f\\x ^ ||9t/||y- As in (i), 

again using smallness of this implies \\g\\x ^ ll^fi'lly- 

(iv) To prove the ^-regularity result, consider the equation dtf + DBf = 0, which 
implies 

n 

\\dj\\y = \\DP H Bf\\ y » W( p -HB)(dif) + Pn{diB)f\\ y 

i=i 

since D = DP-u and the operator D has estimates ||-D/i||2 ~ YH=i ll^^lb f° r & H 
h G D(D) H Ji. (The latter is straightforward to verify with the Fourier transform.) 
Here P-^ denotes orthogonal projection onto H; it commutes with <9j. This yields 
the bound 



\difh « \\(P H B)d t f\\y < \\d t f\\y + \mB)f\\y* < (1 + ||^S||0||/||* < ||/| 



X 



if max(||£9jA||*, ||£d^4||*) < oo, where we used that PuBt '■ H — > H is an isomorphism 
in the first comparison. Conversely, if ||td^4||* is sufficiently small, then 

n 

\\f\\x<\w\\y<£mf\\y+\\tmw\\x), 
i=i 

where the first estimate is by (iii). Using next that Y^7=i 11^*^11* i s small enough, 
this implies \\f\\ x < \\V x f\\ y . 

As in (i) above, these estimates translate to Hc^gHy < \\g\\x an d \\g\\x < ||V x g||y 
respectively. □ 

Proof of Lemma \10.5\ Assume that the coefficients A satisfy || A— A ||* < oo and has 
distribution derivative d t A e L l £(R l + n ; C(C^ +n ^ m )) such that \\td t A\\* < oo. Fix 
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h G X with distribution derivative dth G y. By Theorem 17. 11 \\SAht — S A h t \\ 2 dt — > 
as e — > 0, where 

pt poo 

S e A ht := / r]?(t, s)Ae- {t - s)A E+£ s h s ds - / r/ e "(t, s)A e - {s - t)A E Q -£ s h s ds = 1-11. 
Jo Jt 

Hence it suffices to bound Hc^S^/i)!^, uniformly for e > 0. 

(i) Differentiate I and write 

td t (I) = [\td tV t)Ae-^ A E+E s h s ds- f\t{t-s)A 2 e-^ A E+£ s h s ds 
Jo Jo 

■ f V ?(d s Ae- {t - s)A )K(s£ s h s )ds = [\td tV t + sd sV t)Ae^ A E+£ s h s ds 
Jo Jo 

- [ V ?(t-s)A 2 e-^ A E G + £ s h s ds+ [ r ] ?Ae- it - s)A E+d s (s£ s h s )ds = I 1 -I 2 + I 3 . 
Jo Jo 

Note that in I3 the distribution derivative d s (s£ s h s ) extends its action to test func- 
tions s H- (r)£(t, s)\e~( t ~ s ^ A E Q + ')*(f), for any <f> G %. Theorem 16.51 and Lemma [5751 give 
the estimate 

ll'slb- < \\d t (t£ t h t )\\ y * < (ll^ll* + \\td t £\l)\\h\\ x + \\£\U\d t h\\y- 
To bound we apply Lemma F6. 71 using the bounds 

/t POO 
\(t- s)\ 2 e- {t ~ s)x \sds<t and J \(t - s)X 2 e~ {t ~ s)x \dt < 1, 

which shows H-^Hy* ^ ll^lly* ^ 11^11*11^11^- To estimate I\, we calculate 

(td t + ad.)r,+(t, s) = i^( v y(i^)rj e (t)Ve(s) + V^^iWMVeis) + s Ve (t) v >(s)). 

From this, we verify that \(td t + sd s )r]+\ < XsuppVn? — ^ Hence an estimate as in 
the proof of Theorem 16.51 shows that ||/i||y* < \\£ ||*||^||#- 

(ii) Next we differentiate II and write 



/oo poo 
(td t rj:)Ae-^ A E^£ s h s ds - J t V ; (d s A e -^ A )E^ £ s h s ds 

/OO pOO 
t(d tV ; + d s r ] ;)Ae- {s ~ t)A E r£ s h s ds + / ^\ Ae~^ A E^ sd s (£ s h s )ds 

= Hi + Ih- 

To bound II 2 , we apply Lemma [6~71 using the bounds 

poo ps 

/ \(t/s)\e^ s ~ t)x \sds < t and / \(t/s)\e~ {s ~ t)x \dt < 1, 
Jt Jo 

which shows ||/J 2 ||y» < \\td t £ \\*\\h\\x + \\£ \\oo\\9th\\y. To estimate Hi, we calculate 
t(d t + d s ) V t(t, s) = tr, (±*M(t)rk(s) + vM V ' e (s)). 

The last term is supported on s G (l/(2e), 1/e), t G (e, s — e), where it is bounded by 
et < t/s. Thus estimates as for II 2 apply. The first term is supported on t G (e, 2e), 
s G (£ + e, 1/e) (and another component which can be taken together with the last 
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term) and is bounded by 1. Splitting this remaining term as in ( |3~T]) . it suffices to 
estimate 



X(e,2e)(t)tr]' e (t)e- tA I r] t (s)Ae- sA E^S s h s ds 
o 



< 



2c 



-tA 



y* 

sA f 



/ r) e (s)Ae E 8 s h s ds 
o 




< 



r] e (s)Ae- sA E^£ s h s ds 



< \\£h\\y < \\£ 



using the uniform boundedness of e tA and Lemma 16.91 This completes the proof. 

□ 



11. Miscellaneous remarks and open questions 

(i) The condition iV*(Vt jX w) 6 L 2 implies that Whitney averages ^ JJ w ^ ty ^u 
converge non-tangentially for almost every x, ie. with \y — x\ < at for some 
a < oo, to some u (x) with u belonging to the closure of C^°(R n ) with respect 
t° ||Va;/||2 < oo. Furthermore, t~ l f 1 'V x u s ds converges weakly to V x u in L 2 as 

t — > (compare Theorem I2.3( i)). In particular ||V x ito||2 ^ ||M.(Vt )!C it)||2- This is 
essentially in [331 P- 461-462], where it is done on the unit ball instead of the upper 
half space, and with pointwise values instead of averages, working with u's solving 
a real symmetric equation. However, the result has nothing to do with BVPs, but 
is a result on a function space. 

(ii) Assume that A e L 00 (R/ I ( _ + ' 1 ; £(C (1+n)m )) and that N*(V t ,xU) G L 2 with u 
satisfying (JTJ) in R^ +n distributional sense. Then there exists g E i/ _1 / 2 (R n ; C m ) 
such that 



(44) 



/ / (-AVt^u, W t , x 4>)dtdx 



for all e C °°(R 1+n ; C 



If d UA u(s, x) := (AVt, x u(s, x)) ± for all s > 0, a; G R n , then t 1 J d UA u s ds converges 

weakly to — g in L 2 as t — > 0. In particular ||g|| 2 < ||-/V*(Vt i iw)||2- This is again 
essentially [33] for the unit ball. See [H Lemma 4.3(iii)] for an argument in R^ +n . 
The equality f|4*4l) justifies that g is called the Neumann data. This result has nothing 
to do with accretivity of A, boundedness suffices. Compare again Theorem 12. 3( i). 

(iii) Theorem 12 .4( i) contains a priori estimates on ^-solutions. A natural question 
is to reverse the a priori estimates for such systems. Does a weak solution to 
(CQ) with \\A — Aq\\c < oo and N*(u) G L 2 satisfy ||V$ )!C it||;y < ||i\T*(it)||2? Same 
question replacing N*(u) G L 2 with sup t>0 \\u t \\ 2 < oo. The smallness of \\A — 
Aq\\c, which implies well-posedness of the Dirichlet problem for 3^-solutions, yields 
a posteriori such estimates. It would be interesting to have positive answers a priori 
(i.e. independently of well-posedness) when \\A — Aq\\c < oo. 

(iv) Is there existence of ^-solutions to the Neumann and regularity problems with 
L 2 data under \\A — A \\ c < oo (or even under the stronger J °° UA(t) 2 dt/t < oo, 
where u>A(t) '■= sup 0<s<t \\A S — Ao||oo)? Is there uniqueness under the same constraint 
on A, provided existence holds? Recall that tools such as Green's functions are not 
available here. 
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(v) Same questions for 3^-solutions and the Dirichlet problem with L 2 data. 

(vi) It is likely that 3^-solutions have the a.e. non-tangential convergence property 

for averages: {w ^ y)l JJw(t, y ) u ~* for a ' e - x G Rn and ~* in 

\y — x | < at. This requires an argument which we leave open. 
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